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I.  INTRODUCTION 


Squeeze-f lira  dampers  are  used  in  the  construction  of  modern  aircraft  turbo¬ 
engines  to  allow  its  rotor  to  pass  through  one  or  more  critical  speeds  to 
reach  its  operating  condition  [1].  Typically,  the  damper  is  formed  by  an 
oil  film  which  separates  the  outer  race  of  a  roller  element  bearing  from  the 
bearing  housing.  Due  to  space  limitation,  the  axial  length  of  the  squeeze- 
film  is  in  the  order  of  the  width  of  the  bearing  and  is  effectively  quite  short 
in  comparison  with  its  diameter.  Usually,  the  damper  is  centered  by  a  spring 
structure  which  represents  a  parallel  connection  between  the  bearing  outer- 
race  and  the  housing.  The  lubricating  oil  is  typically  fed  to  either  or  both 
ends  of  the  squeeze-film  through  a  circumferential  distribution  manifold. 

Effective  use  of  the  damper  for  vibration  control  requires  proper  sizing  cf 
the  damper  to  fit  the  flexibility  charac teristics  of  the  rotor  [2J.  Obviously, 
a  minimal  damping  capacity  is  needed  if  the  kinetic  energy  of  the  flexible 
rotor  as  excited  by  a  mass  imbalance  may  be  absorbed  to  avoid  a  cumulatively 
increasing  amplitude  of  vibration  as  the  critical  speed  is  approached.  At 
the  same  time,  if  the  damping  capacity  is  excessive,  for  a  flexible  rotor, 
the  support  bearing  in  effect  becomes  rigidly  attached  to  the  frame.  The 
critical  condition  would  be  shifted  to  a  higher  speed.  The  inherent  damping 
capacity  is  inoperative  because  of  lack  of  motion  of  the  support  point.  The 
vibrational  response  to  rotor  imbalance  would  increase  unchecked  as  the  higher 
critical  speed  is  approached.  Thus,  use  of  damp-mounting  to  control  imbalance 
response  of  a  flexible  rotor  through  critical  speeds  requires  an  optimum 
damping  capacity.  If  the  rotor  ifc  required  to  pass  through  several  critical 
speeds,  the  optimum  damping  capacity  is  a  compromise  among  all  the  modes  c-f 
interest.  To  design  such  a  damper  system  properly,  it  is  necessary  to  carry 
out  a  thorough  rotor  dynamic  analysis  [3]. 

Due  to  inherent  symmetry,  the  response  of  a  centered  damper  to  rotor  imbalance 
is  a  circular  orbit.  According  to  the  lubrication  theory,  there  is  a  kine¬ 
matic  equivalence  between  a  statically  loaded  journal  bearing  and  a  damper 
in  circular  orbit.  Accordingly,  the  theory  of  plain  journal  bearings  can 
be  applied  to  dampers  provided  the  following  rules  are  observed: 


•  The  static  eccentricity  of  the  journal  bearing  is  equivalent 
to  the  orbit  amplitude  of  the  damper, 

•  The  journal  rotational  speed  u>  of  the  journal  bearing  is  to  be 
replaced  by  twice  the  vibrational  rate  2v.  Since  the  vibration¬ 
al  rate  is  also  the  rotational  speed  of  the  shaft  for  imbalance 
response,  one  in  effect  doubles  the  rotational  speed  when 
journal  bearing  data  is  used  to  obtain  damper  characteristics. 

•  The  tangential  force  per  unit  static  eccentricity  of  the 
journal  bearing,  upon  rescaled  for  twice  the  rotational  speed, 
is  equivalent  to  the  damping  force  per  unit  orbit  amplitude 

of  the  damper. 

•  The  radial  force  per  unit  eccentricity  of  the  journal  bearing, 
again  upon  rescaled  for  twice  the  rotational  speed,  is  equiv¬ 
alent  to  the  stiffness  force  per  unit  orbit  amplitude  of  the 
damper,  which  is  additive  to  the  stiffness  of  the  centering 
spring. 

While  the  above-stated  scaling  law  between  a  journal  bearing  and  a  squeeze- 
film  damper  is  basically  sound,  it  is  ironical  that  commonly  accepted  journal 
bearing  data  and  analysis  are  not  suitable  for  use  to  describe  damper  charac¬ 
teristics.  The  commonly  used  Ocvirk-Gurabel  analysis  for  short  journal  bearings 
(4,5)  is  not  a  suitable  model  for  the  squeeze-film  damper.  This  is  because 
a  journal  bearing  usually  receives. its  make-up  lubricant  supply  through  a 
feed-hole,  which  is  within  the  bearing  surface;  whereas  the  lubricant  flow 
in  a  damper  must  come  in  from  either  end.  In  order  to  remedy  the  situation, 
Shaker  Research  Corporation  proposed  to  undertake  a  study  to  advance  the 
technology  of  squeeze-film  dampers  by  establishing  the  appropriate  analysis 
for  its  characteristics  in  a  manner  consistent  with  the  typical  configurations 
and  lubricant  supply  method.  The  work  was  sponsored  by  U.S.  Army  Research 
Office  under  Contract  Number  DAAG29-78-C-0027  for  the  period  from  July  15, 

1978  to  July  15,  1981,  and  was  performed  under  the  guidance  of  the  Technical 
Representative,  Dr.  Edward  A.  Saibel  of  the  Engineering  Sciences  Division, 

U.S.  Army  Research  Office. 
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II,  SUMMARY  OF  MAJOR  FINDINGS 

The  theoretical  study  of  a  squeeze-film  damper  undergoing  a  circular  orbital 
motion  was  performed  in  a  manner  which  is  also  applicable  to  a  statically 
loaded  journal  bearing,  so  long  as  the  lubricant  supply  is  fed  through  either 
end.  It  was  established  that  the  film  rupture  process  must  satisfy  the  flow 
continuity  condition  throughout  the  "cavitation  domain" .  The  Swif t-Stieber 
condition  [6,7]  would  govern  the  breakup  boundary,  whereas  the  fill-back 
boundary  is  established  according  to  the  flow  balance  between  surface  trans¬ 
port  (through  the  cavitated  domain)  and  the  re-established  film.  It  was 
further  concluded  that  the  short-bearing  approximation  originally  advanced 
by  Ocvirk  and  Dubois  [4]  can  be  adapted  to  satisfy  the  above-stated  con¬ 
ditions  around  the  film  rupture  boundary.  Closed  form  solutions  were  found, 
from  which  a  full  set  of  numerical  results  was  compiled  and  presented  in 
the  dimensionless  form.  The  following  general  conclusions  can  be  surmised. 

•  The  state  of  a  ruptured  film  is  defined  by  three  parameters; 
namely,  the  pressure  in  the  ruptured  domain  (or  the  largest 
orbit  size  without  rupture),  the  operating  orbit  size,  and 
the  axial  location  of  the  rupture  axis. 

•  The  pressure  in  the  ruptured  domain  is  an  unknown  parameter 
of  the  system.  It  can  be  fixed  at  the  ambient  level  only  if 
a  "drain  away  condition"  is  maintained  at  the  discharge  side. 

•  If  lubricant  is  allowed  to  drain  away  at  both  ends  of  the 
damper,  the  lubricant  depletion  process  occurs  with  the  evo¬ 
lution  of  a  succession  of  "necking  down"  patterns. 

•  If  the  damper  is  driven  into  an  orbit  of  an  ever  increasing 
orbit  amplitude  after  film  rupture  has  commenced,  the  overall 
extent  of  rupture  would  increase  accordingly  and  can  reach  into 
the  convergent  portion  of  the  gap.  As  a  consequence  of  the 
enlarged  rupture  domain,  the  damper  would  be  more  "spring-like". 

•  With  a  free  draining  discharge  at  one  end,  an  appropriate 
supply  pressure  fed  to  the  other  end  would  effectively  fix 
the  operating  characteristics  of  the  damper. 


* 
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Experimental  observations  were  made  by  direct  viewing  of  the  rupture  pattern 
through  a  transparent  damper  bushing.  It  was  found  that  minor  details  in 
the  outboard  flow  passageways  can  prevent  free  draining  of  lubricant  dis¬ 
charge,  When  this  happens,  the  rupture  process  appeared  to  be  associated 
with  local  nucleation  sites  and  developed  slowly  with  the  aggregation  of 
void  pockets  which  fluctuate  as  the  dynamic  load  passes  by. 

The  unsteady  motion  of  a  rupture  boundary  was  also  treated  theoretically. 

This  was  done  in  a  manner  consistent  with  the  lubrication  theory.  Flow  in 
the  ruptured  domain  was  modeled  as  an  adhered  film.  Inertia  and  body  forces 
were  neglected  as  well  as  meniscus  surface  tension  and  interfacia]  phase 
transfer.  A  constant  "cavity  pressure"  was  assumed  but  was  allowed  to  be 
either  higher  or  lower  than  the  film  pressure  bordering  the  rupture  boundary. 
A  "law  of  independent  action"  was  proposed  so  that  the  thickness  of  the 
adhered  film  and  the  sweeping  rate  of  the  rupture  boundary  may  be  separate 
interfacial  parameters.  The  commonly  used  Swif t-S tieber  condition  was 
shown  to  be  a  special  case  of  the  "law  of  independent"  action  for  the  qui¬ 
escent  time- independent  state.  Computation  of  the  adhered  film  field  in 
the  ruptured  domain  would  have  to  proceed  simultaneously  with  the  computation 
of  the  film  pressure.  The  former  was  governed  by  a  real  characteristic  in 
the  world  diagram  of  the  sliding  motion  and  would  have  to  be  reconciled  in 
the  determination  of  the  sweeping  rate  of  the  rupture  boundary.  The  new 
theory  can  be  used  conjunctively  with  any  combination  of  the  surface  sliding 
velocities  and  with  the  presence  of  the  squeeze-film  action.  It  can  be  used 
is  the  basis  for  the  numerical  simulation  of  large  amplitude,  motion  of  a 
ruptured  film,  with  or  without  sliding,  including  the  non-circular  orbit  of 
the  squeeze  film. 
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An  Improved  Short  Bearing  Analysis 
for  the  Submerged  Operation  of 
Plain  Journal  Bearings  and 
Squeeze-Film  Dampers1 

H\  dlltui  ing  the  Jtlm  pressure  to  assume  some  subambient  value  and  by  alloti  .tut  natural 
boundaries  of  the  film  to  form  in  the  unloaded  region,  the  short-bearing  tht  ttr\  of  tH'firh 
and  Duhois  is  extended  to  m  elude  11  detail*  d  desert  fit  inn  of  the  em  it  atom  zone  i't.  >>  ait* 
natu cavitation  eon  figurations  ar e  slum  n  to  be  possible .  rendering  different  recent  n»  (/\ 
(orbit  size ).  a'titude  angle  (phase},  for  the  same  load  and  minimum  film  pressitn  T},e 
first  configuration  features  an  enclosed  cavity  maintained  at  a  subambient  tend  and  is 
called  "O"cai"tation.  u  hich  is  crudely  emulated  b\  the  conventional  half-film  "approxi¬ 
mation.  The  second  configuration  features  ambient  level  side  cavities,  the  boutidanc'  of 
which  are  draun  inside  the  bearing  ends  by  a  sub-cavity  film  pressure,  and  is  call *  d 
cavitation  The  "/  cavitation,  which  is  initial*  (/  b\  the  aggregation  of  entrained  bt.bf  b's 
in  the  ambient  fluid ,  can  present  itself  in  the  form  of  multiple  striatums  causing  substan¬ 
tial  loss  of  load  capacity 

I  Introduction 

The  short  I  tearing  approximation  is  a  powerful  analytic  al  device 
for  the  estimation  of  journal  bearing  forces.  Keren!  interest  in  the  use 
of  squeeze* film  clamper  for  vibration  control  has  drawn  considerable 
attention  toils  use.  The  geometrical  constraint  of  a  small  (!./!)•  rat  io 
turns  out  to  be  commonly  quite  acceptable.  The  conventional  short  - 
bearing  lormul.ie  lor  the  Itearing  forces,  however,  are  also  based  on 
an  additional  ht  uristic  hypot  bests  that  the  subambient  port  ion  of  tie 
film  pressure  Mould  not  l>e  realized.  Under  the  disguise  of  this  •* rea¬ 
sonable”  mocel  concerning  film  pressure,  a  gross  inconsistency  con¬ 
cerning  the  lubricant  circulation  pattern  is  often  overlooked.  By  ne¬ 
glecting  the*  vthamhienf  film  pressure,  axial  flow  at  the  l>ea ring  ends 
can  only  lie  net  ward.  In  actual  practice,  the  flow  leakages  through  the 
bearing  ends  ire  usually  mad-?  up  hv  the  supply  of  lubricant  through 
a  feed  hole  in  the  unloaded  portion  of  the  (tearing  gfip.  Howeve  r,  the 
journal  (tearing  may  also  be  used  in  a  submerged  arrangement  The 
squeeze  film  damper  often  operates  in  this  manner.  In  a  submerged 
arrangement,  lubricant  end  leakage  from  the  loaded  f  tort  ion  of  the 
bearing  film  must  lie  made  up  by  the  re-entry  of  the  lubricant  into  the 
bearing  gap  on  the  unloaded  side.  Thus,  the  presence  of  some  mi- 
bambient  pressure  in  the  bearing  film  cannot  Ite  avoided. 


1  Supported  *»y  the  U.  S.  Armv  Research  Office  under  Contract  Number 
1)AA(i»-7M'-miJ7,  Project  Number  P-1A657-K. 

Contributed  hv  ihe  Lubrication  Division  of  THE  AMERICAN  Society  OF 
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Kver  since  the  formulation  of  the  lubrication  theory  by  Osltorne 
Reynolds  |]|,  the  tpiestion  of  suhambient  film  pressure  has  been  a 
matter  of  much  controversy.  (Uimhrl  |2|  pion<*ered  the  popular  rot  ion 
that  an  adequate  approximation  is  attained  by  discarding  that  portion 
of  ihe  full  solution  of  the  film  pressure,  which  is  subambient.  Afu  r 
the  short  bearing  approximation  was  introduced  |3).  it  Itecnme  im¬ 
mediately  linked  to  (UimM's  idea,  the  short  bearing  film  is  usually 
taken  to  l>e  the  convergent  portion  only.  Meanwhile,  more  precise 
treatment  of  the  film  cavitation  process,  which  is  presumably  ihe 
mechanism  for  preventing  the  film  pressure  to  fall  much  Mow  the 
ambient  level,  had  l>een  attempted  in  various  ways.  Swift  |4|  and 
Shelter  |5|  inde|>endenlly  revived  Reynolds’  idea  that  at  a  film  rupture 
boundary.  the  gradient  of  the  film  pressure  must  vanish.  Jako  tsson 
and  Klolterg  |6|  demon.si rated  the  procedure  to  treat  the  rupture 
boundary  in  a  two  dimensional  field.  Ktsion  and  Pinkus  |7|  ga\e 
computed  results  including  a  detailed  acoont  for  the  feeding  slot  or 
hole.  Wakurai  et  al.  |8|  considered  a  journal  bearing  with  a  circum¬ 
ferential  feeding  slot;  they  introduced  the  idea  that  a  short  bearing 
type  analysis  can  be  combined  with  the  treatment  of  the  film  rupture 
a’ndition  of  Swift  and  Stieber.  While  Jakohsson  and  Ploherg  br  >ught 
out  the  significance  of  a  suhambient  condition  through  an  assumed 
vaftor  pressure,  the  recent  efforts  |7,  H|  still  dealt  with  situations  where 
the  film  pressure  could  not  he  suhamhient.  The  cavitation  bou.id.tr> 
conditions  were  also  re-examined  with  respect  to  local  flow  fields, 
including  the  roles  of  the  separation  streamline  |9|  and  the  surface 
tension  in  the  meniscus  |I0|.  Floherg  returned  to  the  question  of 
minimum  film  pressure,  this  time  via  the  concept  of  liquid  tensile 
atrength  |ll|,  he  sought  to  determine  the  pattern  of  str  atum 
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streamers,  wh  ch  have  lieen  Men  from  photographs  of  the  cavitation 
domain.  In  ctnirast  to  the  rupture  condition,  Kloherg's  stri.it ion 
streamer  analysis  depends  on  the  assumption  of  zero  net  flow  across 
the  boundary.  Nau  examined  experimental  evidence  of  cavitation  on 
face  seals  [12)  He  concluded  ihat  there  is  a  distinction  lietween  the 
"separation  ci  vi tat  ion”  which  is  associated  with  a  direct  communi¬ 
cation  between  the  cavitation  domain  and  the  ambient,  and  the 
"ligament  cavitation"  whit  h  features  an  enclosed  cavitation  domain. 
Interest  in  vibration  control  by  means  of  squeeze  film  dampers  led 
to  experiment  il  information  (tertaining  to  the  submerged  operation 
of  short  squeeze  journal  lieariitgs  |12,  U|.  Both  White  and  Tonnesen 
found  abrupt  changes  in  the  hearing  characteristics.  Marsh  |  lf»|.  re 
ferring  to  White  s  data,  suggested  that  sulidivision  of  the  fluid  film 
can  take  place  to  cause  n-ductton  of  the  dynamic  load  capacity  of  the 
damper. 

It  is  clear  thai  Cumber*  simple  model  is  not  always  adequate  to 
cope  with  cavitation  phenomena  in  industrial  machine  elements.  It 
also  became  evident  that  one  may  not  he  able  todescrilie  all  cavitation 
phenomena  in  fluid  films  with  a  single  model  such  as  the  Swift  -St  ieber 
lor  Reynolds)  boundary  condition,  even  if  one  is  willing  to  remain 
ignorant  alxnit  processes  which  jiertain  to  local  flow  fields.  At  the  mine 
lime,  although  toe  short  hearing  formulae  as  derived  with  (iumliefs 
hypothesis  ar*  not  suitable  for  a  submerged  journal  liearing.  there  is 
no  fundamental  reason  to  preclude  the  use  of  the  short  liearing  ap¬ 
proximation  with  a  more  appropriate  treatment  of  the  cavitation 
process.  The  work  of  Wakurai.  et  al  already  demonstrated  the  fea 
sibility  of  treating  an  anthi«  nt  tavitation  with  the  Reynold  *s  condition. 
The  present  work  represents  another  step  in  this  direction. 


The  respective  description  of  the  film  thickness  is.  for  the  >  atu 
displacement  of  a  rotating  journal 

h  *  (’  —  e  cos  0  fta» 


and  for  the  circular  orbit  motion  ol  a  squeeze  film  damjier 

h  =  ( '  —  c  cos  (W  —  i*t )  liiht 

Definitions  of  the  dimensionless variables  for  the  two  circumstances 
are  indicated  l»elow: 
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The  film  pressure  is  measured  above  the  ambient  reference  In  tin 
case  of  the  damper  in  a  circular  orbit,  the  equivalent  sliding  vi  ctoi 
is  in  the  opposite  direction  with  the  amplitude  2rR.  thus  resulting  t 
the  corresponding  definitions  of  /*  and  rj. 

rfhe  short  liearing  approximation  presumes  Of  Ok  iH  1  Ol*' Ok  i  and 
d/d»|  (H * d/Vdij)  to  be  of  a  similar  magnitude  so  t hat  in  the  limit  A  * 
0.  equation  (21  i*  reduced  to 


0k\  on  on 

For  the  axially  uniform  film.  OH  0{ 
lui  ion 


II  =  \-V>  i  I 


=  0.  equation  1  1 l  lias  th  ■  ■ 


II  Review  of  Background 
Convention]  Short  Bewaring  Analysis.  With  the  approximation 
of  a  constant  viscosity,  the  laminar,  incompressible  lubrication  theory 
applicable  to  a  thin  fluid  film  in  a  cylindrical  anuulur  gap  is  repre 
sented  by 

1  0  I  ,0/A  0  {u)Rh  hA  0p\  Oh 

- \h '  —  +  - — - (1 ) 

12m  >;  \  tel  2  \'2»R  Of* I  01 

For  two  circumstances,  (a)  the  static  displacement  of  a  rotating 
journal,  imd  t‘n  the  steady  circular  orbit  motion  of  a  squeeze  film 
damper,  Equation  (l)  can  he  reduced  to  a  single  dimensionless 
form: 
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[\lu  11 1 )  are  invariant  with  respect  to  j*  and  are  determined  from 
suitable  boundary  conditions.  The  most  commonly  accepted 
boundary  conditions  are  associated  with  the  imposition  of  the  am  >ient 
pressure  of  the  ends:  i.e. 


\u< «  ±  n  *« 


so  that 

11 1  *  0;  ii(>  *  — ~ 
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—Nomenclature—— 

An  *  consUmi  of  integration  in  the  neckdown 
region  for '  cavitation 
C  *  radial  charance 
I)  *  diametei  of  journal 
v  m  eccentric  t\  of  journal 
F  *  function  of  m  given  by  equation  127) 
F#,  f’  T  m  dim.-iisionless  radial  and  lungenlial 
force  comp  »  ents  equation  (17) 

//  *  dimensionless  luhruanl  film  thickness 
or  liearing  jap,  h/C 

h  *  lubricant  film  thickness  or  liearing  gap 
L  •  liearing  t  xial  length' 

I *  *  dimensionless  film  pressure;  pi'V 
ifipuiR2)  f<  r  static  inurnal  hearing,  p("2/ 
( \ZpvR'2\  for  squeeze  film  damper 
p  *  film  pressure  above  ambient 


I  dll 


R  *  journal  ladius,  l)f  * 
t  *  time 

?  *  axial  coo  H  nate 


«»„  *  fratiion.il  fouette  trans|H>rt  across 
neck -down  boundary 
(  *  operating  eccentricity  ratio,  c/r 
tf  -  eccentricity  ratio  for  cavitation  inci¬ 
pience 

f  *  dimensionless  axial  coordinate.  2 zfl. 
fu  -  dimensionless  axial  distance  from  film 
center  to  cu\  ital  ion  boundary 
ij  *  generalized  circumferential  coordinate 
/)  -  circumferential  angular  coordinate 
X  *=  slenderness  parameter,  /.//> 

M  *  viscosity  otefficient  of  fluid 
i*  ■  frequency  of  oscillation  of  damper 
II  ■  short  liearing  dimensionless  pressure. 
Fil)/t.)2 

Ilo,  1 1 1  ”  integration  constants  in  generalized 
short  lieuring  solution 
Hr  *  value  of  II  in  cavity 
Ilmm  *  lower  Itound  of  II  for  fluid  film 
X  *  Sommerfeld  angle. 

.  L  /ITT  fa)? 
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Subscripts  for  if 

U  =  w  idest  point  of  "tf*  cavitation 
b  *  break  up  Isaindary  in  "0"  cavitation 
bt  *  break-up  initiation  in  "tf*  cavitation 
<  i  »  cavitation  incipience  (Miint 
/  *  fill  Iwck  liotindary  in  "0"  cav  tatiou 
\  *  minimum  pressure  point  in  7**c«  vita 
turn 

u  *  neck -down  Iniundary  in  "/"  cavitation 
m  *  neck  down  initiation  in  "I"  cavitation 
n  *  widening  boundary  in  *7"  cavitation 


Subscripts  for  H  and  Q 

etc.  *  Hi  iff,),  Qiiff,)vU\ 


Subscripts  for  Fp  and  Ft 

FF  •  full  film 
(  !)  ■  cavitation  domain 
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and  thus,  the  so-called  “short  tearing  solution"  is  given  hv 


II 


mm 


£  !I(  <  0 


i VM 


n--i(i-r-iQ  »«» 

This  equation  shows  that  the  lubricant  film  pressure  would  te  su- 
hamhient  in  the  interval  0  <  i)  <  x.  Since  the  latter  cnndit  ion  is  gen¬ 
erally  believe  *d  to  te  unlikely  in  an  indust  rial  environment,  tin  “short 
bearing  solution"  is  often  accompanied  by  the  “half-film  approxi¬ 
mation"  by  ietting 

Il(i)  <  »j  <  jr)  *  0  (91 


The  “half  Him  approximation"  removes  the  |»erfect  antisymmetry 
with  respect  to  rj  of  II.  Accordingly,  the  global  characteristics  of  the 
lubricant  as  defined  by 


Fr  ~  ^  dC  ^  II  cos  rj  dij 

Fr  *  -  J  df  £  llsinifdif 
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are  drastically  altered.  In  the  case  of  the  "half-film  approximation." 
the  above  integrals  are  calculated  only  in  the  interval  -w  <  »?  <  0.  The 
integrated  results  are  well  known: 
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Rationale  for  the  Half  Film  Approximation.  Fixation  of  l he  Him 
inlet  at  q  =  —  7r  and  the  exit  at  i}  *  0  is  a  fair  approximation  for  most 
industrial  situations.  The  exit  location  is  justified  on  the  assumption 
that  the  lub 'leant  film  cannot  sustain  any  significant  suhainhient 
condition.  With  the  lubricant  fe«*dhole  placed  in  the  unloaded  jHirtion 
of  the  bearing  gap.  the  inlet  film  would  be  formed  not  too  far  from  the 
maximum  gt>p  While  it  is  realized  that  the  inlet  location  would  vary 
somewhat  dep*  nding  oil  the  specific  feeding  arrangement,  the  overall 
film  force  does  not  depend  on  the  precise  inlet  location. 

Sum  of  end  leakages  and  the  exit  flow  are  made  up  by  the  flow 
through  th  inlet  film.  Any  excess  in  the  flow  supply  through  the 
feedhole  would  te  discharged  to  the  amhient  in  the  unloaded  gap. 
Thus  the  totil  flow  through  the  feedhole  ultimately  reaches  t he  am¬ 
bient  in  the  form  of  axial  flow  at  the  tearing  ends. 

Peculiarity  of  the  Submerged  Bearing.  In  a  submerged  tearing, 
because  then  it  no  separate  feedhole,  the  net  flow  through  the  tearing 
ends  has  to  I  e  zero.  The  end  flow  in  the  loaded  region  is  necessarily 
outward,  which  must  te  balanced  by  inward  end  flow  in  the  unloaded 
region.  The  |  recise  flow  balance  of  flow  determines  the  extent  of  the 
load  carrying  *11ra.  Some  degree  of  subambient  condition  is  un¬ 
avoidable.  1  hi  desired  operating  characteristics  of  a  submerged 
tearing  may  very  well  require  the  lubricant  film  to  develop  a  signifi¬ 
cant  subambient  level.  Regardless  of  the  magnitude  of  the  suhamhient 
condition,  the  structure  of  the  fluid  film  on  the  unloaded  side  will  have 
a  strong  influence  on  the  overall  film  force.  If  some  form  of  cavitation 
should  develop,  the  fluid  film  in  the  unloaded  gap  should  te  narrower 
than  the  tearing  width.  Thus  one  should  expect  the  short  tearing 
approximation  to  te  at  least  equally  applicable  in  the  unloaded  gap 
as  in  the  loaced  gap. 


Ill  A  New  Model  of  Cavitation  for  a  Submerged 

Bearing 

The  conceptual  discrepancies  in  the  conventional  short  tearing 
solution  (“half  lilm  approximation"!  can  te  overcome  by  the  following 
hypotheses. 

1.  The  lubricant  film  can  sustain  some  Htihamhient  and  even 
tensile  condition  but  is  txamded  by  a  lower  limit 

II  a  1  { m«n;  Ilm.nCO  (121 

2.  The  film  pressure  in  the  cavitation  region  is  governed  by  the 

range 


Implied  along  with  the  second  condition,  inertia  and  capillarv  ef¬ 
fects  would  te  neglected.  The  general  short  tearing  solution,  ne 
giecting  misalignment,  as  given  hv  equation  (5!,  will  be  retained 

One  must  also  give  some  attention  to  the  conditions  at  the  -  av  tia 
tion  boundary.  The  commonly  accepted  Reynolds'  condition  tor  the 
lubricant  to  enter  the  cavitation  boundary  is  (grad  II  **  0  at  II  *  lift. 
While  no  explicit  description  regarding  the  flow  Held  inside  the  cav 
itation  region  is  associated  with  the  Reynolds'  condition,  it  is  usiiatlv 
presumed  that  the  velocity  profile  is  a  straight  line  in  appro  idling 
the  I  mi  ind  ary,  but,  upon  crossing  over,  lubricant  is  transported 
through  the  cavitation  region  in  a  thin  layer  which  adheres  lo  the 
moving  surface.  Since  pressure  gradient  is  negligible  within  the  cav¬ 
itation  region,  the  lubricant  cannot  te  redistributed  in  the  direction 
transverse  to  the  surface  motion.  To  satisfy  continuity,  the  film 
thickness  at  the  re-establishment  boundary  (relative  to  that  at  the 
break  up  boundary!  must  be  consistent  with  the  pressure  gradient 
which  would  require  a  concavity  in  the  velocity  profile.  This  de¬ 
script  ion  is  illustrated  in  Fig.  1. 

Strictly  s|ieaking.  according  to  such  reasoning,  the  cavit.itcd  domain 
is  not  capable  to  transfer  shear  stress  from  one  surface  lo  the  other. 
However,  it  is  a  common  practice  to  invoke  an  alternate  mole!,  in 
which  fluid  transport  through  the  cavitated  region  is  contained  m  a 
multitudeNif  st nation  streamers.  Thus  the  cavitation  boundary  ac¬ 
tually  is  not  smooth,  but  displays  a  fluctuating  pattern  correspi aiding 
to  the  spacing  of  the  striatum  streamers.  The  fluid  contained  in  the 
st  nation  streamers  develops  little  pressure  gradient,  but  contributes 
to  a  transfer  of  shear  tetween  the  tearing  surfaces.  Figure  2  shows  an 
illustration  of  cavitation  streamers. 

One  may  further  indulge  in  the  dichotomous  description  of  the  fluid 
structure  (in  the  event  of  cavitation!  in  the  global  scale.  In  the  case 
of  a  submerged  plain  journal  under  load,  the  condition  of  equation 
( Id!  can  te  satisfied  in  two  distinct  configurations  The  central  region 
may  break  up  and  give  way  to  an  isoharir  rpgion  at  11/  *  IIm  n  <  0. 
The  central  cavitation  region,  however,  is  serrated  from  the  ambient 
ends  by  “side  bands,"  through  which  transverse  pressure  drop  from 
the  ambient  continually  feeds  the  lubricant  with  make  up  flow  Such 
a  configuration  for  cavitation  will  te  referred  to  as  the  *  0"  petti  rn. 
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Fig.  i  rim  cavitation  of  Hw  “O”  c  iwHgiif  IHw 


One  may  retain  the  ''dichotomy"  of  the  flow  structure  at  the  mi*  ro 
scale.  That  is.  Reynolds'  condition  for  film  break  up  can  be  imposed 
where  applies bU ,  while  the  si  reamer  concept  can  be  retained  to  render 
a  conservative  estimate  for  the  friction  torque.  Alternately,  side 
boundaries  of  the  unloaded  imrtion  of  the  fluid  film  may  recede  from 
the  bearing  end*.  Presumably  gaseous  constituents  will  fill  up  the  two 
side  regions  between  the  bearing  ends  nnd  the  necked  dow  n  film.  The 
pressure  in  the  necked  down  region  will  be  subambient,  but  the  decree 
of  necking  down  ensures  the  condition  llmin<KUllr*0.  The  la’ ter 
form  of  cavitati  m  will  be  designated  “/.**  The  lower  bound  value  of 
pressure  may  require  the  film  to  be  subdivided  into  a  multiple 
pattern  consisting  of  more  than  one  film  strip.  For  the  present  pur 
p«ee,  symmetry  about  the  mid -plane  prevails  in  either  case,  therefore 
it  is  sufficient  to  limit  our  attention  to  one  half  the  axial  extent,  i.e. 
0  <  i.  In  the  case  of  a  multiple  “/**  pattern,  $  *  0  would  he  the 
centerline  of  each  film  strip.  Figs  it  and  4  show  the  film  patterns  of  the 
two  configurations;  the  respective  analyses  for  the  two  cavitations  will 
be  dealt  with  separately  in  subsequent  sections. 

IV  Analyses  of  the  “0”  Cavitation  Configuration 

Cavitation  ca  i  commence  at  the  midplane  when  the  lubricant  tilm 
enters  into  the  divergent  gap  s<  that  the  midplane  pressure  reaches 
the  condition 

Qh. 

IK?*.  C  -  0)  «  -  y"  I*'  <0  M4> 

where,  tfh  local**  the  initiation  of  the  "break-up"  lioundary.  and 

<*.  «<*<**> 

_ «  sin  ifa 

1 1  -  t  ro*  ):* 

Clearly,  given  II/*,  this  condition  can  exist  only  if  t  exceeds  a  lower 
bound,  which  can  he  called  the  cavitation  incipience  eccentricity  ratio 


Fig.  4  T‘  cavitation  configurations,  single  and  wndtiple  atriationa 


Fig.  9  Eccentricity  ratio  and  cantor  for  cavftation  Indplanca 


(for  the  "0"  pattern!,  rj ll<- >.  When  f,  is  being  exceeded,  cavitation 
begins  as  a  spot  which  centers  around  the  location  of  the  minimum 
pressure  of  the  full  film  solution.  The  location  of  the  spot  from  wh  ch 
cavitation  grows  can  be  called  the  cavitation  incipience  center.  tri. 
which  is  the  root  of  the  stationary  point  of  (1  in  the  divergent  port  on 
of  the  full  film  solution.  U,.  »?, « l  as  functions  of  Il(  are  given  by 


(,  Sill  !).  , 


*  Il<  *.  t,  cos  *  -  |\  1  +  24i,2  -  l\  tiro 
2(1  —  t,  cos  Tft  t  )•*  4 


They  are  shown  as  graphs  in  Fig.  5. 

Suppose  that  r  >  e,  for  a  given  Ilf,  the  film  break-up  location  ip* 
can  be  solved  from  equation  (141.  then,  for  ip.  >  ip*.  the  break  up 
boundary  )  is  to  be  found.  The  short  bearing  solution.  Kq  nation 
(51  is  applicable  for  fn  $  C  ^  I.  The  boundary  conditions  are 
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ip.  is  the  circumferential  coordinate  of  the  break-up  boundary. 
Furthermore.  Reynolds*  cavitation  condition  requires 
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and 
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Making  use  *1  equation  •  i4i. 


1  -  v;-- 

(h 


i  lilt 


<2m 


The  stationers  cundil  urn  of  Cu  turns  out  to  coincide  with  the  l.*ca 
tu»n  cd  the  minimum  pressure  »f  the  full  film  solution,  i.e. 


ilif  fu 


\\  I  4  24 1-  -  ti 


c21 » 


q,t  locate*  the  v  ulrst  point  of  the  "O"  cavity  or  the  narrowest  point  o| 
the  "side  bam:,  which  terminates  the  "break  up”  boundary.  ren 
dermg 


(«(,,  m  4<>Ufu>  IteprrseiHadw  will  tains  of  the  break  up  region  ot  the  *0” 

cavitation  are  •  hown  in  Kig  6 

For  ^  1  hi*  "ode  Iwitd*  widens  and  equation  (17)  no  longer 

applies  Instead  the  sum  if  the  transverse  pressure  flow  ind  (he 
transport  thus  ghthe"<»’’e  ivitv  should  become  the  required  lull  film 
fow.  which,  in  the  c a  >1  ih-  *)u>rt  liearing  theory,  is  simply  (he 
Couette  flow  os  illustrated  in  Fig,  7,  the  necessary  condition  is 


n  p-  1 1  -  «  cos  If/ 1 1  Ol!’ 

iitf,  t  ci*  rjt,  -  cos  H(\  IdtT/'o 


The  circumlerenhat  o**rdmate  is  designated  m  to  indicate  the 
"fill-buck”  prison.  Together  w  ith  the  pressure  boundary  condit  ions 
equation  4  Ihl.  the  short  liearir  g  theory  now  reduces  to 


( Ucos  fjj, 


t|f>“  f  sin  TfjU 
dr)f 


c«>2 

*  “  7(1  —  #  COS  If/I^llr 


1*241 


which  describe*  (lie  geometry*  i  f  tl»e "fill  bat  k”  boundary.  The  upjMU 
limit  of  r\f  is  re  idled  when  Co  teduces  to  zero  again. 

Kquation  42  U  «.  an  or  dinar*  differential  equation  and  can  be  ,n 
tegrated  for  ij/  ?  to  troin  the  initial  condition 


C<>4  Vf  *  »jii)  *  tfolmai  t2M 

which  is  given  by  equation  4221.  ifb  is  a  function  of  <fo  according  to 
equation  420).  therefore,  the  ditferential  equation  is  rumlinmr,  and 


Fly.  •  Break-up  boundary  ef  tba  **•"  aaeWadan  canftguratkM 


its  rigorous  solution  can  only  lie  obtained  numerically.  Howevei.  a> 
shown  in  Fig  ti,  lor  most  part  of  Co  <  4 solm«».  fib  deviates  \ery  b  tl« 
from  >ih,  Therefore,  one  would  esjied  that  linearization  ot  equation 
1241.  obtained  by  replacing  n*  with  a  suitable  constant  rtf,.,.  car  Ih* 
quite  satisfactory*  A  reasonable  choice  for  this  constant  is  given  b> 

1  /•'tnW 

C«IS  * - I  COS  >2ti) 

which  define*  »he  correct  total  "four He  transjKirt"  across  the  entire 
bre.»  ktu*  boundary  The  In  e  iri/.ed  problem  has  a  closed  form  snlu 

tioii 


4  1  -  u.l 


where 


H„- H,.. 

n  -  nK 


|  1  ~  t  («■  ht,4t  j 


21I( 


t/f  -H^) 


Ih'(rtf)  -  F0\  1 27 1 


Ft  ip  I  -  1 1  4  -  i  *  jiff  -  ;t«  1 1  4  j  I**|  sin  rtf 


a  ,  i  , 

+  -r  sin  2  *t, - r*  sin  t 

4  12 


An  =  Ft  It,  *  I till 

Ho  *l*i  cos  t) a;  H s.,  *!—i  cos  r tb.. 


The  almve  relations  Idr  the  "fill  back”  domain  holds  for  t  u.)ni ,,  * 
»*..(  ni  t  *  0  l  b*  latter  i  iMldit  nm  is  reached  at  rj,  *  rj, .  where  a  full  I II in 
,u  ross  t  he  entire  hearing  w  nil  h  is  re-established.  r)T  is  solved  nuttier 
icatly  from  equation  4  27)  bv  setting  Co  *  t). 

Thin  approximate  solution  can  be  called  the  "mean  transport  ap 
proximal  ion."  The  full  cavitation  pattern  c  alculated  according  to  thi- 
approximation  is  shown  for  various  combinations  of  («,.  r)  in  Fig.  K 


V  Analysis  of  the  **  P  Cavitation  Configuration 

The  new  model  for  cavitation  proposed  in  Section  111  admits  ,t 
second  cavitation  configuration.  Suppose  the  fluid  film  remains  intac  t 
at  the  midplane  while  the  film  pressure  is  maintained  above  the  lower 
Uamri  Nmin  by  drawing  in  liquid  gas  interfacial  menisouses  from  the 
ends.  Such  a  ph  ’  condition  is  quite  reasonable  if  gaseous  mb 
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I 


Flf.  7  CuMiwdly  Mrow  a  IW-bacfc  boundary  tw  a ’t"  captation  cawWfwra- 
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V 


F*.8  Boundary  pattern  ol  “0"  cavitation  configuration  Fig.  §  Continuity  octom  th#  nocliing-down  boundary  m  on  r  cavitation 

configuration 


stance  is  allowed  to  cumulate  near  the  edges  of  the  hearing  film,  The 
corresponding  short  hearing  solution  (neglecting  surface  tension ) 
is 


I]  =  y  <sV-  i'!)Q„  (28) 

applicable  in  a  yet  uns|»ecified  range  0  <  ?|  <  jr.  Jot  rjn)  is  the  libation 
of  the  men  st  us  tmundary.  and  Qn  m  Q(ijn),  *)«  being  the  circumfer¬ 
ential  coorJi  late  of  the  necking-down  boundary.  It  is  important  to 
recognize  th.  t 


(29) 


varies  linearly  with  the  distance  from  the  midptane  and  assumes  its 
largest  amplitude  at  the  meniscus  Naindarv  for  anv  value  of  jj,,.  This 
edge  film  pressure  gradient  is  responsible  for  feeding  the  “necked- 
down”  film  which  centers  at  the  midplane.  A  crucial  question  to  be 
settled  at  t  hit  point  is  the  degree -of- flu  id  transport  in  the  side  voids 
once  the  'Tec king -down”  process  has  started.  This  issue  does  not  exist 
with  the  ”0”  configuration  because  the  Reynolds’  cavitation  condition 
was  impost  d  m>  that  dl  I  /dj  vanishes  at  the  break-up  boundary.  Since 
axial  fluid  ransport  in  the  side  voids  cannot  be  sustained,  the  overall 
circumferential  fluid  transport,  including  both  that  through  the 
“necked -down”  film  and  that  through  the  side  voids,  must  be  in¬ 
variant  once  the  necking-down  process  starts.  At  a  differeni  ml  >-cale 
the  flow  b*  lance  condition  across  the  necking-down  boundary  is  de¬ 
picted  in  F  g.  9.  Accordingly,  the  axial  coordinate  of  the  necking  down 
boundary  must  satisfy 


(I  -«„>//„ 


HR  I 


Substituting  into  equation  (28)  and  setting  {  =  0,  the  center  lire 
pressure  becomes 

lit f  ■  o.  if, i  *  -  -  A*2  — -  j Hn * 

2  d»|„  / 

The  right-hand  side  has  a  minimal  condition  given  iiy 

4«2  cos*  rj ,v  +  i  cos  —  5c1'  *  0  <  .(4  i 


where  rjv  designates  the  location  of  the  minimum  film  pressurt 
Solving  for  i  cos  rj*.  one  finds 

1  _ _ 

(  cos  rj\  *  -  |\  l  +  80<  -  -  1|  I :t  “» I 

8 

One  can  then  calculate,  writing  (^\  for  Qins I 


lit iT  *  0,  tf\)  *  < :w> • 

One  may  olwerve  that 

*I\<B..:  Us<Q>,  *  IT. 

Consequently,  incipience  of  necking-down  depends  an 

•  ■  nun 


represents  ‘'fractional  Couette  transport”  at  the  necking  down  Mr 
boundary.  Since  the  side  voids  represent  the  interruption  of  ambient 

feeding  of  the  subambien'  film,  they  would  be  dried  up.  Therefore  1 

the  correc  value  is  nr,  *  0.  Accordingly,  the  continuity  condition  which  is  essentially  the  same  as  the  incipience  of  breaking-up  in  th** 
across  the  necking  down  laaindary  is  '*0”  cavitation  process.  An  important  characteristic  of  the  **1”  envi 

j|\  tation  is  that  the  extent  of  necking-down  is  finite,  as  sikh.  a*  th** 

(“J  *  ”  H„  j-  (.10)  condition  for  its  initiation  is  satisfied.  In  contrast,  upon  im  ipiencc 

'  of"0"  cavitation,  the  extent  of  cavitation  grows  grad  ua  IK  liefinnin*: 

Making  u*e  of  Kquai  ion  (29)  and  recalling  the  definition  ol  (^.  one  with  an  infinitesimal  amount.  The  neck  down  initiation  p  tint  v< 
obtains  the  differential  equation  for  •  coincides  with  ?>,,  for  c  *  r  The  corresponding  centerline  pressure 
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is  the  lower  Itoun  l  value.  A><  exceeds  7,,  moves  toward  the  mini 

mum  gapwhih  i-ie  i  or  respond  ing  centerline  pressure  stavs  at  II . 

until  i)f|l  reacln  >  7\  7,,,  -  q\  is  determined  h>  the  <  midition 


,  *  sin  7,,, 

Quit,.  I  *  ~  ~  -  Hu 

( l  -  «  Cits  7„,  I 


l.Wl 


with  turther  in  rcase  ot «.  7,,,  t  ills  Itelow  »j\  while  the  minimum  film 
pressure  stays  u  7\ .  so  i  hot 


7  v  i  =  V  - 


Us 


Together  wi  h  equation  *d2  .  tme  can  calculate 

tHs\ 

'.•'ll/.* :i  r..<i»s*  i— j 


Since  Hn  ca  unit  i*e  smaller  than  th<*  minimum  gap  H\  7  :z  (t  =  1 
—  f.  it  follows  that 


S'a  7,  )  < 


1  - 


0»<7O 


(41  i 


If  the  right!  and  side  of  equation  (41  )  is  larger  than  unity.  hen  the 
reck  down  ini'  iat  ion  p  <int  7„  >  0  can  lie  I  mind  as 


7„  *  cos*1  [l  -  H.x  si»l f?.v >]/e| 


<42 


for  the  solut ion  >  f  a  simple  "l"  ct ^figuration.  If  t he  right  hand  >ich 
of  equation  <  41  )  should  he  less  then  unity,  then  the  maximum  width 
of  the  film 


_  Hsttdnsi 
I-  < 


<42 


Widening  ol  the  tilm  t hi n commences.  Because  the side-v  oids are 
dry,  a  relation  s:  nilar  to  equa! ion  (82)  holds,  i.e. 


// .  mi'  7u  1  *  const.  *  ,4M 


•Fr.  Ft)  '■ 


(WH,  WT) 


D 


m/.\  2 

7. 


=  2  f  Icon  e,  “  sin  f  IldC 

Jw  •/« 


hoih  from  the  standpoint  ot  computational  convenience  ami  I  mu 
the  standpoint  of  physical  interpretations.  diMerent  procedures  to- 
the  computation  of  him  force  components  are  dev  t*ed. 

“IT  Configuration.  Fur  the  “0**  configuration,  it  is  usef  il  to  miIi 
divide  the  film  force conifioncurs  into  those  in  the  “full -fi!m”d  at  .in 
covert'd  by  7,  <  tj  ^  ip,,  and  those  in  the cavilated  domain  co\<  -e<:  h. 
7*.  *•:  7  €  tfr  +  2ir.  Thus 

( Fit ,  F  /  )  =  t  f  /,■  F I  Iff  4-  1 b  /, .  F  I  h  t ,  Is 

Closed  form  expressions  are  found  lor 
•  f'H  hi 


t  l«M 


<  (cos  rji,,  —  cos  rpitcos  rj,.  +  cos  ijr 


-  t  cos  T)h,  COS  7-  ^ 


ih'rhy 


where. 


M  1  —  t  cos  i/f,  )-(  l  -  i  cos  t),  l  • 


<  <  1  -  (  i  ;  I  1 

— — — ^ - -  t\f,.  -  y.i-  -tsm2\>„  -siii7\.u  *  an 


{s  only  a  fraction  of  the  h*  arin  :  length  The  spare  length  would  then 
!h*  occupied  b\  one  or  more  additional  film  stripts)  as  illustrated  in 
Fig.  4 

The  neck inp-( low n  process  terminates  at  the  maximal  condition 
of  //„,  which  occurs  at  7,,  =  x  yielding 


|WN\ 

w»(  7'.  31  t-i<  *)  * - hitijv)  till 

\l  +  il 


<  4>* 


vu  is  the  oircu  nlercntial  coordinate  of  the  widening  houndary  ami 
//„  *  Hiiju  1.  It  is  clear,  to  maintain  continuity.  p,( tju  =  7rl  **  ^,,<7. 
!S  x)  and  H„  ti axl-//,;(iy,  »  ir). one  must  require  A u  =  A„  He 
tause 

//„  *  Hni'>*  -  17, 4  1  (4<h 

the  geometry  of  y0(ij„  )  is  an  exact  mirror  image  of  pit  7,,  t  about  the 
line  of  centers  1'he  domain  of  the  widening  process  is  n  *  rj,.  5?  Jjr  - 

I’m  • 

An  essential  rhstirctive  leat  ire  which  is  common  to  all  simple  uul 
multiple  *7”  >  «.  ifigurat ions  s  the  precise  symmetry  of  ihe  film 
iHiundaries  <am!  am  isymtnetrv  of  the  pressure  field)  with  respt  ct  to 
the  line  of  cent*  rs.  I'lus  is  however  in  part  due  to  the  omission  ot 
curvature  in  tl  e  1  i r«  uniterenli.il  vehn  itv  profile  by  the  short  licaring 
theory. 

VI  Film  Force  CnlcuUtion 

I'he  film  force  can  )w  calculated  in  terms  of  its  radial  and  tangi  ntial 
(omponents  relative  t<  the  line  of  centers  In  accordance  with  the 
iummiin  geonu  rioil  eon  vent  inn,  they  can  l>e  expressed  nondi 
mensionallv  a> 


,  1  L  •  l  +  I  7 

\  -  2  t  m  1  \  - Ian - 

V  1  *  '  2, 

In  the  cavitatcd  domain,  one  may  use  II,  =  II  as  the  referi  1 
pressure.  From  equation  I  IK) 

j”  .11  -  II,  ht< 

=  J'1  hi  -  ii,  ur 

k\»r  the  break-up  region,  equation  1 19)  applies,  theretore 
/• 1  I 

|  (II  -  ll«  ui\  -  -  *  1 1  —  Collie  i°r  ns,  *  70 

In  (he  till-back  region.  p,( 7,  t  is  given  by  equation  t*27t.  Thus. 

b'i  )fi, 

*  2  |(sin  7r  “  sin  7 ),  (cos  —  cos  7^  )|I1< 

2ll,-  /*** 

- ; —  I  (1  -  s'o  I  (cos  tih,  -  sin  Uhlaiih 

d 

u  -  s>*  Q 


J%  tr  ♦  » 

1 1  -  .Co*  + 

*K< 


«  Hr  j 

x  loos  7;,  —  sin  fyh/fy  .V!.» 

which  cun  l>e  readily  computer!  by  numerical  quadrature  with  tht  aid 
of  etpiat ions  (19)  and  127). 

**l”  Configuration.  Consider,  first,  the  simple  "I"  configuration. 
For  the  “full-film'*  domain,  - 7„,  <  7  ^n(.  Kquat ions  (49  and  inui 
apply,  upon  replacing  (7*,.  7. 1  hv  I  if*,,.  -  7n  h  resj>ectivel>  .  Tht  s 

=  0 


iF  1  hr  s  ;  <11  -  *  ’Cr-’  |v„  -  -  sin  2\ni|  »-"*** 


where. 


■  2  tan  1  V  - tan  — 

r  1  -f  2 


For  the  necked  down  region,  7,,,  ^11,^1,  front  equation  128). 

•  iio 


J*  1  f  iia 

lids'-  I  llc/C 

1)  «/o 


(47) 


Cm1 

*  -  - —  t  .V- » 

Because  of  the  antisymmrtrical  character  of  the  film  pressuie.  t 
is  not  necessary  to  treat  the  widening  region,  r  <  ijh.  <  2x  -  7^.  *r*p 
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tern  Fig.  10  Geometrical  parameters  of  ”0”  cavitation 


VII  Result!} 

The  incip  ence  of  cavitation  is  governed  by  the  lowest  pressure 
which  can  he  -uMaiicd  by  the  lubricant  film.  It  is  convenient  l*» 
translate  this  threshold  condition  into  the  cavitation  incipience  ec¬ 
centricity.  (,  as  given  in  Fig  5. 

When  the  o|»erating  eccentricity  exceeds  the  incipience  condition, 
(  >  r,,  the  film  >n  the  unloaded  side  can  break  up  to  form  a  ca\  ilation 
pattern  of  th-*  "0”  configuration.  The  essential  geometrical  parameters 
of  this  pattei  n  are  the  I«h -at ions  of  the  breakup  and  re-establishim  nt 
points  and  the  width  friction  of  the  cavitation  domain.  They  are 
functions  of  <  >  r  I  as  shown  in  Fig.  10.  The  break-up  point  is  always 
located  in  the  divergent  part  of  the  (rearing  gap.  A>  <  >  (i  increases, 
the  break -u )  |K»int  accordingly  shifts  toward  the  minimum  gap; 
meanwhile.  1  a*-  re  establishment  point  extends  into  the  large  gap  area 
and  the  width  of  the  cavitation  domain  assumes  a  larger  fraction  of 
the  tearing  I  *ngth.  The  re-establishment  |>oint  ultimately  would  pass 
the  maximum  gap  and  w«hjI<I  he  located  in  the  convergent  pan  of  the 
liearing  gap  The  "far  extending  cavitat ion  condition."  wit  b  l he  film 
re-establishi  H  in  the  convergent  gap,  can  occur  even  at  a  very  "modest 
operating  eccentricity  if  f,  is  small,  for  instance,  at  I*  -  (VI.  r,  * 
0.1165). 

The  axial  ir**ssure  profile  for  the  "0"  cavitation  configuration  tan 
tie  classified  into  three  main  gn>u|is  as  illustrated  in  Figure  1 1  In  the 
full  film  regi  m,  ij,  £  n  *  ija,.  the*  familiar  symmetrical  paraliolu  pro 
files  apply.  I  he  i  enfer-line  pressure  is  above  ambient  in  the  con¬ 
verging  film  and  is  negative  in  the  divergent  film.  In  the  break-up 
region.  7a»  5  7  ^  7<».  t hi  parabolic  profile  is  split  by  the  cavitation 
domain,  and  it  has  a  zero  slope  at  the  break-up  boundary  according 
to  the  Reyw  Ms'  boundary  c< mdition  In  the  fill  hac  k  region  while*  the 
ends  of  the  pressure  profile  arc  fixed  respectively  at  the* »  av  ly  and 
ambient  pre  sure>.  il*  shape  may  lie  conc  ave  or  conve  x  depending  on 
whether  it  is  located  in  the  divergent  or  in  the  convergent  |h>H<oi>  of 
the  hearing  ran.  In  fact.  for  the  combination  of  a  large «  and  i  small 
it  is  possible  for  the  pressure  to  la*  above  ambient  so  that  the  film 
may  lie  losing  tluid  t*»  teth  'he  cavity  and  the  ambient. 

The  overall  film  lores-s  t  »r  the  "O'*  cavitation  configuration  are 
shown  grapl  ically  in  Figs.  1 1  and  1*1.  In  Fig.  12.  the  tangential  furrr 
for  various  \alues  of  t,  are  plotted  against  ♦  together  with  the*  "coin- 


[ 

i 


DIVCI 


COHVtatCHT  0AP 


Fig.  1 1  Axial  protura  prcWti  tor  "0"  cavitation 


plete  film**  ami  the  "hall  film"  results  for  comparison  For  all »  *.  • 
the  "complete  film"  curve  is  approached  as  cavitation  is  either 
suppressed  or  is  of  an  insignificant  extent.  On  the  other  hand,  tor 
moderate  values  off,  iff,  is  less  than  0.0f»,  the  tangential  force  an  te 
significantly  smaller  than  the  value  hased  on  the  "half  filn"  ap 
proximal  ion.  Thus,  depending  on  the  value  of  the  minimum  t  Im 
pressure,  which  is  associated  with  the  value*  off,  according  in  Mgurr 
5.  the  "half  film"  tangential  force  may  Ik*  in  error  by  n  file  lor  «t  two 
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Fig.  12  fore*  for  “0"  cavitation 


either  way.  In  Fig.  lit,  the  radial  force  is  compared  with  the  "fc.df  f  *- 
result.  Aga  n.  for  <  <  r,.  the  radial  force  along  with  th«  tfliu 
cavitation  Incomes  negligible.  For  very  small  values  of ;  *•  *•  -  •. 
force  can  h  larger  than  the  “half  film"  result.  In  short  ■'* 

solution  may  be  deficient  either  by  overlooking  the  ca<  + 
pience  threshold  or  by  disregarding  the  “far  extend  in  s  :  ».■  a 

condition.’ 

The  “ / ”  •avitation  proci  HS  may  take  place  inste  >r»  the  *  >at 
tern,  with  t  be  film  lioundaries  receding  front  the  ring  ends?  o  o  rm 

a  necked -d  >wn  centra’  dim.  Depending  on  the  magnitude  of  the  op 
erating  em  ntricity  t  relative  to  that  for  cavitation  in«*ipie»i«e.  a  variety 
of  film  patterns  may  develop.  These  possibilities  are  de  <•*  bed  below  . 
As  t  first  exceeds  film  neck-down  commences  at  if.,,  *  rj<;  according 
to  equation  1 15)  where  the  center-line  pressure  is  lowest.  The  film 
continues  10  neck  down  as  the  gap  increases  with  the  film  width 
varying  inversely  as  Ihc  gcp.  The  neck  of  the  film  is  located  at  the 
maximum  ;*npt  lieyond  which  the  film  width  increase's.  The  film  ge 
ometry  is  s\  mine  tried  about  the  line  of  centers  while  the  tilm  pressure 
is  antisymnelrical,  consequently  the  film  force  is  purely  tangential. 
Because  su  dia  necking  do  vn  process  causes  an  abrupt  change  in  the 
film  pattern  as  (  exceeds  f,,  the  film  fierce  accordingly  decreases 
abruptly.  As  '  l>ecomes  larger.  .  moves  toward  the  minimum  gap. 
carrying  w  th  it  the  minimum  center  line  pressure  which  is  at  the 
lower  boun  I  value  As  if(IJ  p  isses  below  rj/v,  which  is  given  by  equation 
C15),  with  furl  her  increase  efr,  the  minimum  center  line  pressure  t  still 
maintains  at  the  lower  hound  value)  stays  at  17 >;  the  center  line 
pressure  at  i]  „  is  higher. 

Ultimate  Is  (he  neck  down  initiation  point  would  reach  the  mini¬ 
mum  gap.  I’ii  M  her  inert  ase  of  4  requires  the  appearance  of  a  narrower 
secondary  ^tiip  which  <  on’ acts  the  main  film  strip  at  the  minimum 
gap.  Minimum  pressure  p  tints  are  fixed  at  tf\  in  both  strips.  The 
minimum  pressure  in  the  main  film  strip  stays  at  the  lower  hound 
value  whih  1  lat  in  the  sec  mdarv  strip  is  higher.  With  increasing  t. 
the  widths  of  the  main  and  secondary  strips  respectively  decreises 
and  increaM  '  until  they  become  equal.  A  still  larger#  would  then  re 
quire  anothi  r  secondary  strip  to  appear,  starting  a  new  cvcle.  While 
the  film  pattern  takes  on  such  variations  as  illustrated  in  Fig.  1.  the 
axial  pressure  profiles  remain  as  simple  parabolas  with  the  only 
stipulation  that  the  minimum  pressure  remains  no  less  than  the  lower 
hound  value  'Hie  film  lore,  which  is  purely  tangential  for  the  entire 
family  of  the  patterns,  is  shown  in  Fig.  14.  For  <  <  *,,  the  I  ibricant 
film  is  complete,  the  corre ^ponding  film  force  is  shown  as  t  dashed 
line.  As  1  it  aches  r,.  the  film  force  ahrupily  drops  to  corres|*md  to  the 
abrupt  change  of  the  film  pattern,  the  neck-down  region  new  spans 
rjr,  $  rj  4  2  it  -  ij(l.  Further  increase  of  t  curiously  is  accompanied  liv 
a  reductiot  of  the  film  ton  e.  This  trend  continues  through  ihe  dmi- 
hle-striatio  1  condition.  When  the  double -atriat  ion  pattern  approaches 
the  condition  of  strips  of  equal  width,  the  film  force  a<  (nelly  rises 
somewhat  with  *  hut  turns  downward  again  after  the  third  strip  has 


u. 


Fig.  13  Radial  tore*  tor  ”0“  cavitation 


Fig.  14  FHm  tore*  tor  **f"  cavitation 


lH»en  formed.  Again  it  resumes  the  upward  trend  when  all  dm  e strip- 
approach  an  equal  width,  then  down -turn  again  with  the  formation 
of  the  fourth  stria!  ion  strip.  This  priKTss  continues  indefinite  y  u  :  il- 
4  increases  toward  unity.  The  numeral  at  the  right  margin  dedguMe 
the  number  of  striation  strips  in  the  film  pattern.  One  mav  observ.- 
that  for  a  given  the  film  force  is  always  smaller  than  tha  ai  ihi 
cavitation  initiation  condition  (r  *  in  the  case  of  the  “/“conligu 
rat:nn.  For  the  “0“  configuration,  in  contrast,  the  film  forces  ahv.iy  • 
increase  with  «  despite  the  growth  of  the  cavitation  domain  as  shown 
in  Figs.  12  and  hi. 

VIII  Dtacuaaion* 

*  The  present  work  extended  the  applicability  of  the  short  tearint 
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approximation  to  treat  cavitat  on  domain*.  Irrespective  «if  the  short 
hearing  approximation.  however.  important  questions  art*  raised  re 
garding  the  exi>  tome  of  suhamhient  film  pressure  and  the  appropriate 
choice  of  cavil. ilion  iHiund.iry  conditions. 

Allowance  01  a  suhamhient  pressure  is  essential  in  the  analysis  of 
submerged  journal  lH*anng  in  order  that  the  overall  continuity  of  lu 
hricant  flow  can  he  satisfied,  Iakohsson  and  Flnherg  had  correct  I v 
staled  the  cast*  |4>(  and  presented  results  which  are  equivalent  to  the 
“0"  cavitation  problem  Floh.*rg  further  recognized  that  the  How 
structure  in  th  ■  cavitation  domain  is  not  smooth,  hut  comprises  of  a 
number  of  striMion  streamers;  he  associated  the  number  of  stru  tion 
streamers  will  the  tensile  strength  o)  the  lubricant  jllj  Flolierg. 
however,  restricted  the  latter  )>oint  of  view*  to  a  scale  which  is  much 
smaller  than  tl  e  «verall  cavitation  domain.  Nau  distinguished  sepa 
ration  cavitation  from  ligament  cavitation  |12|,  which  roughly  cor 
r/spond  to  the  *7Mand  "0"  captation  configurations  treated  in  the 
present  work.  The  present  point  of  view  may  be  regarded  as  a  con 
solidation  of  th  .■x*  ideas,  The  p  edibility  of  a  suhamhient  pressurt  and 
the  proper  cho  ce  ol  cavitaf  iun  Ixiundary  condit  ions  cannot  be  sepa 
rated,  l’he  tot  d  pin ure  consists  of  (he  collection  oi  the  following 
points: 

1.  The  cavitation  Imundery  can  l>e  identified  as  a  break  up 
boundary  or  a  ie  establishment  lx>undary  according  to  it s  geometrical 
relationship  wifh  the  sliding  vector.  In  crossing  a  break  up  boundary, 
the  sliding  vect  *r  leases  .»  flaid  film;  while  crossing  a  re  establishment 
boundary,  it  enti  rs  a  fluid  film. 

2.  In  approaching  a  break-up  boundary,  if  the  film  pressure  i> 
higher  than  tht  pressure  in  the  cavitation  domain  then  it  must  a* tain 
a  zero  gradien  i.|xm  ri  aching  the  breakup  Itoiuularv  This  cor  re 
sponds  to  the  Swift  Sli*d»cr  or  the  Reynolds* condition. 

3.  If  the  li  n>  pressure  is  Mibcavity  in  the  neighborhood  ol  the 
break-up  boundary,  a  mat '.ere  pressure  gradient  must  exist  at  the 
boundary  such  that  t  here  is  no  net  flow  across  it.  An  exception  ft  ilii'- 
rale  is  when  th<  break-up  hnin  lary  has  become  parallel  to  the  sliding 
vector  and  is  turning  intoji  re  establishment  boundary. 

4.  Fluid  in*  m  port  across  the  cavitation  domain  follow's  t h«*  sliding 
vector  and  miet  lie  conserved  lietween  the  break-up  and  the  re  es 
tiblishment  p  ints.  Ifllie  break-up  boundary  borders  on  a  sultcnvitv 
film,  then  the  cavitation  domain  must  lie  a  dry  void  and  the  r»*  es 
tablishment  hmirdary  must  he  followed  by  a  film  pressure  rise:  again 
there  can  be  no  ret  flow  across  such  a  re -establishment  boundary. 

Note  that  the  above  conditions  concern  only  thin  film  parome 
ters. 

Inaccuracy  a ssociated  with  the  use  of  the  short  bearing  approxi 
mation  are  present  in  the  internal  field  and  along  the  cavitation 
boundaries.  The  field  inaccuracy  is  associated  with  the  omissiun  ot 
the  circumferential  Poiseudle  ’omponent of  the  film  flow  relative to 
the  Couette  component  Its  magnitude  can  readily  be  estimated  bv 
computing  and  comparing  these  components  from  the  solution  of  the 
film  pressure.  In  erms  of  flow  cate,  the  field  inaccurac y  in  o|/.//>{  but 
in  terms  of  film  pressur*,  The  field  inaccuracy  alsouf  ects 

the  break-up  boundary  of  "O’  cavitation  directly  in  that  there  »s  a 
residual  circumfi  rential  pressure  gradient,  which  amounts  to  o|/./7>; 
when  scaled  h\  t  ne  axial  pi essure  gradient  At  the  re-establishment 
boundary,  bee  uise  the  •  in  umferential  pressure  gradient  would  fa 
cilitate  filling  up  th «•  gap.  t he  I  oundary  would  shift  more  toward  lh« 
maximum  gap  than  the  Tvrt  bearing  result.  From  this  point  of  view 
one  might  expet  t h  it  the  defu  lency  of  the  "half  film"  approximation 
is  prolxably  ho  newliat  overestimated  in  the  present  work.  Foi 
cavitation,  the  circumterentifd  pressure  gradient  would  upset  th< 
precise  symmetry  between  the  necking -down  and  the  widening 
boundaries.  The  slo|>e  of  the  necking-down  boundary  should  be  re 
duced  while  t  int  of  the  widening  Ixuindary  should  l»e  increased 
(Consequently,  the  pressure  field  would  not  be  precisely  anti-sym 
metrical,  so  th  it  (he  film  fore*  would  have  a  radial  component. 

A  more  disconcerting  featuie  of  the  short  bearing  approximation 
is  revealed  in  tht  pressure  fiel  l  of  the  film  with  "O'*  cavitation.  The 
situation  is  charly  illustrated  in  Fig.  11.  l'he  suhamhient  pressure 
prttfile  change,  from  tht  simp  e  paraltola  in  the  full  film  to  the  split 
parabola  in  the  break  up  region,  'l'he  split  parabolas  still  have  a  zero 


slofie  at  (lie  cavitation  b  darv,  ami  (he  pressure  profile  negoti.itt  - 
the  transition  smooth**,  •  pm  entering  the  fill  back  ngioii.  il 
pressure  profile  at  tains  a  finite  slope  at  the  txnmdury,  still  the  change 
takes  place  smoothly.  This  happy  situation  unfortimaleh  is  no  ,i- 
sured  everywhere.  The  cavitation  fkniiidarv  keeps  the  film  pressure 
at  the  lower  tiound  value.  Where  the  re-establishment  Ixumdarto  join 
each  other  at  the  center  line,  the  following  pressure  profile  resit nu  >. 
the  full  film  character  of  a  paraltola  with  zero  slop*  at  the  tenter  hi 
in  contrast  to  the  finite  slope  at  the  neighboring  re-estahlishnent 
boundary.  The  situation  is  actually  a  little  worse.  The  center  lii.e 
pressure  of  the  re-established  full  film  is  always  higher  than  the  lower 
Umml  value;  in  fact,  it  can  Ik-  above  ambient  if  t  he  re-establishiuent 
point  is  in  the  converging  portion  of  the  gap.  Thus,  not  only  sum  *th 
ness  but  also  continuity  of  the  pressure  profile  is  \  iolated  in  I  hi  n r 
cumferential  direction  at  this  particular  location.  Such  an  anoi.uilv 
is  usually  associated  with  the  solution  of  a  singular  perturbation 
analysis  |It>|  which  reduces  the  i»rder  of  the  differential  field  equ:  ’ 
when  the  asymptotic  limit  of  an  infinitesimal  numerical  paramet*  r 
is  carried  out.  l'he  short  (tearing  approximation  is  an  applicate  n  ol 
the  singular  perturbation  analysis.  T'he  solution  is  an  asi  ruptoti,  jb\ 
valid  solution  everywhere  except  on  the  line  of  film  re-esiablishn  »-ttt 
In  this  case,  the  deficiency  is  the  jump  in  the  pressure  profile.  V 
“boundary  layer."  type  correction  to  the  short  Ikm ring  solution  can 
lie  found  by  the  method  of  matched  asymptotic  expansion  (I6|.  The 
circumferential  extent  of  the  “lioundarv  layer"  is  0|/,//)|.  Wit  hit  th- 
“boundary  layer."  a  differential  liquation  of  the  full  order  applies.  I  hc 
“boundary  la\er"  correction  smoothly.  all>eii,  rapidly  blen<ls  the  tw.« 
pressure  pr**files  found  in  short  hearing  approximation  on  either  sine 
of  the  re -establishment  |tuml.  One  may  note  t  hat  Itecause  the  w  ultii 
of  the  “boundary  layer"  is  0|/.  /i|.  the  deficiency  in  the  film  t  irir 
found  by  the  short  Itearing  approximation  is  also  limited  u>0{/.  l)\ 

The  two  cavitation  configurations  are  both  mathematically  \  ;ilid 
solutions  (irrespective  of  (he  short  Itearing  approximation!.  F.xperi 
ments  of  White  and  Tom  use  n  on  squeeze  film  damjters  1 13.  14}  m 
defiefidefit/v  lend  some  credence  to  the  possibility  of  two  di/iereol 
modes  of  operat ion.  It  is  therefore  useful  to  query  what  ph\sica;cir 
cumferences  may  pr«»mote  the  presence  of  one  over  the  o*her?  H  • 
due  to  the  answer  is  found  in  I  he  concept  of  dry  voids  which  elii.g  t<* 
the  sides  of  the  necked-down  (ilmtsl  in  the  configuration,  ('he 
(Hissihle  presence  of  these  voids  would  he  assm’iated  with  the  alum 
dance  of  entrained  gas  in  the  supply  fluid  and  a  stagnant  slate  ii  the 
immediate  neighborhood  of  the  bearing  ends.  The  squeeze  film 
damper  is  more  likely  to  operate  with  the  “/*'  cavitation  configure!  *>n 
than  the  rotating  journal  I >eca use  journal  rotation  is  likel*  to  cause 
enough  local  circulation  to  prevent  the  formation  of  the  voids.  In  am 
case  a  high  through  flow  rate  without  excessive  agitation  it<  avow!  ga- 
entrainment  >  is  probably  a  feasible  approach  to  forestall  thp 
configuration  which  would  cause  a  substantia!  reduction  in  -  he  film 
force. 

The  natural  short  Ixearmg  cavitation  parameter  is 


pw'fr 

T  he  scaling  parameter  (/.-/>(.  h<mever,  may  not  lie  involved  if'  th.* 
physical  priKvss.  Therefore  it  is  useful  to  examine  the  results  w  th  a 
fixed  value  of  the  “universal  cavitation  parameter" 


Mwff- 

insiead  of  llmin.  One  may  also  take  the  opportunity  to  further  scru¬ 
tinize  the  significance  of  the  short  hearing  approximation  h>  »  »tn 
paring  the  results  for  “0"  cavitation  with  the  results  of  .Iakohsson  and 
Floherg  |6|.  The  curves  in  Figs  15  and  16  show  such  comparison^  for 
two  values  of  the  universal  cavitation  parameter 
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in  terms  of  the  eccentricity  Incus  and  the  load  capacity,  ren>ecti'-elv. 
As  compared  to  the  half-film  approximation,  present  analysis  brings 
out  the  lack  of  cavitation  at  small  eccentricities  as  indicated  h\  th*> 
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Fig.  IS  Efftds  o4  UD  and  cavitation  parameter  on  eccentricity  locus  In  “0” 
cavHatton 

flat  portion  * >1  the  locus  near  the  origin.  This  discrepancy  of  the  half 
film  approximation  is  more  prevalent  for  a  small  IJI)  The  results  of 
.Jakohsson  nnrl  Flolterg  |ti|  are  in  very  good  agreement  with  the  present 
results  so  f;  r  is  the  errentr icily  locii  are  concerned.  For  the  smaller 
cavitation  parameter,  although  the  locii  still  leave  the  origin  in  the 
horizontal  direction,  they  rapidly  turn  around  and  fall  under  the 
half-film  remits  The  load  capacity  curves  in  Figure  lt>  reveals  the 
deficiency  <  f  the  hall  film  approximation  in  estimating  load  <  apa«  itv 
throughout  the  small  and  moderate  values  of  eccentricity,  *  <  <Mv 
Present  results,  expressed  in  terms  of  the  second  power  law  (with 
respect  to  /  //>).  shows  that  the  effect  due  to  the  dependence  oft  on 
!JD  is  moderately  prominent  in  the  range  0.2  <  t  <  O.fi  for  the  ‘wo 
values  iL*D  --  0.2*>,  U.M  considered  at  the  larger  cavitation  parameter 
(-0.1  Oh  Tl  e  results  of  Jal.obsson  and  Floherg  indicate  a  stronger 
effect  of  UD  At  L/D  =  0.2;-.  the  short  lairing  approx i mat  ion  is  quite 
satisfactory .  At  L/D  =  O.ft,  the  short  hearing  approximation  overes 
ti  mates  the  load  input  itv  bv  12  percent  up  toe  -  0  fv  at  r  =  O.tihy  2ti 
percent,  an  I  at  t  »  n.A  by  .»'*  percent.  At  the  smaller  cavitation  pa 
rameter  (-t'.oi ),  the  load  capacity  is  significantly  reduced  by  the  '  far 
extending’*  cavity. 

IX  Con<  fusions 

The  present  analysis  has  shown  that  the  level  of  subambient 
pressure  ie\  el.  which  can  hi  sustained  hv  the  fluid  film,  is  an  impor¬ 
tant  parameter  which  controls  the  operating  characteristics  of  a 
submerged  ournal  hear  ng.  The  commonly  used  "half  filin'*  solution 
turns  out  to  hr  a  crude  approximation  for  the  solution  w  ith  ca citation 
in  I  he  “0**  <  mfigurat  «>r. 

With  "if  c  avitation,  t h*  lower  hound  of  film  pressure  sets  the 
pressure  lew!  in  hi  taitnted  domain  ind  determines  the  .uno  int 
of  fluid  feei  ing  fr  *m  the  ambient  to  the  fluid  film.  With  a  low  ca\  if  v 
pressure  (large  *,  i.  there  is  a  copius  amount  of  feed  flow  and  the  full 
film  condition  can  he  read  Iv  re-established,  keeping  the  extern  of 
cavitation  to  a  small  area.  \N  ith  the  cavity  pressure  only  slightly  below 
the  am  hie  iv  (-mill  r  ».  a  far  extending  cavitation  region  would  exist, 
the  full  film  re  establishment  point  mav  reach  deeply  into  the  con¬ 
verging  por  iriiofthc  bearing  gap.  In  bet  ween  these  extreme  condi¬ 
tions,  the  "luff  film"  solution  represents  a  fair  approximation. 

"/”  caviu  f  •  »n  causes  the  minimum  pressure  to  exist  at  an  internal 
point  of  the  fluid  film.  A  requisite  condition  for  the  existence  of  this 
cavitation  configuration  is  toe  aggregation  of  gas  or  vapor- filled  voids 
in  theambiriv  fluid.  Th  ■  fluid  film  IsHindaries  recede  away  from  the 
lienring  end  ^  n  I  lowing  dry  Voids  to  iNinpv  the  wide  gap  portion  of  the 
end  regions  of  the  hearing.  I'he  void  extent  may  lw*  stretched  to  the 


Fig.  IF  EfUcIt  of  UD  and  cavitation  pressure  on  load  capacity  in  0 
cavitation 

full  circumference,  thus  preventing  feeding  of  fluid  film  from  th< 
ambient  entirely.  In  addition.  the  voids  may  in  tors  per  si  narrow  I  i  In. 
strips  to  form  a  multiple  stiintion  pattern.  The  occurrence  ot  -  !' 
cavitation  is  accompanied  bv  a  significant  reduction  of  the  fit  id  f'lm 
force  and  thus  may  be  regarded  as  a  pathological  mode  of  >p*  ration 
With  the  short  hearing  approximation,  the  fluid  film  force  ‘or  /' 
cavitation  is  purely  tangential. 

A  generalized  cavitation  model  allows  the  two  different  ca\itntii*n 
configurations  to  exist.  The  new  model  is  a  consolidate *n  <»»'  the 
boundary  condition  of  Swift -St ic!ht  for  film  rupture  and  hat  <*t 
Flofierg  for  striatum  streamers.  Hither  one  of  these  conditioi  s  i*  .d 
lowed  to  prevail  depending  on  the  film  pressure  level  mar  th< 
boundary  relative  to  that  in  the  cavitation  domain 

Successful  useof  the  short  Injuring  analysis  in  conjoin 'ion  wit )  iht 
two  cavitation  configuration-  underscores  the  inqMirtaut  pro  r<  quoin 
of  satisfying  overall  flow  continuity  of  the  lubricant  filmprio  t>  tin 
determination  of  the  film  prosure  Held  and  the  associated  fluid  him 
force.  Imposition  of  the  short  bearing  approximation  is  nec«s»uri!\ 
accompanied  by  the  presence  of  abrupt  lines  in  the  pressure  firWi~  'or 
the  "0”  cavitation  pattern  since  t  he  short  I  rearing  analysis  is  mh  -nut  Iv 
a  singular  perturbation  technique.  The  approach  presented  lien 
however,  is  much  more  realistic  than  the  oonventinnn'  "hal  film’ 
short  bearing  approximation  and  jiermiis  further  improvement  by 
the  technique  of  matched  asymptotic  expansions. 

Present  results  for  “O'*  cavitation  agree  closely  w  ith  earlier  re-ults 
of  Jakohsson  and  Klolrerg  tor  /.//)  =  0.2ft  in  the  entire  range  of  ec 
centricitv.  For  L/D  -  0.,».  the  agreement  is  still  within  2d  pener'  for 
all  <  <  <Uv 

Hrfcmtcw 
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.discussion. 


I.  Etsion3 

The  pa|>er  tie.Us  .in  important  problem  which,  as  was  pointed  out 
correctly,  has  been  a  matter  ol  much  controversy  ever  since  the  lor 
mulation  of  the  lubrication  theory.  The  submerged  Itearing  is  in  many 
aspects  similar  to  a  face  seal.  Indeed  the  concept  of  the  "0”  cavitation 
configuration  i>  known  in  the  seal  literature  since  the  work  by  Findlay 
|17|  in  1968.  Findlay  usi-d  the  -ame  approach  as  in  the  present  work 
to  show  a  cav  it  t  ton  region  <  xtemling  into  the  converging  film  of  a  face 
seal.  The  rationale  liehind  tht*>e  analyses  is  the  required  net  zero  end 
flow  in  the  case  of  submerged  operation.  The  ”0”  cavitation  concept 
preserves  this  requirement  and  maintains  flow  continuity  on  the 
cavitation  bou  alarms  as  well.  It  results,  however,  in  a  pressure  dis 
continuity  across  the  fluid  film  at  tj  =  ryr.  Another  possible  solution 
which  preserves  both  the  net  zero  end  flow  and  the  pressure  continuity 
can  be  obtained  from  equation  (81  by  setting  11  *  II  at  f  »  As  This 
solution,  however,  does  not  c«  mply  with  the  flow  continuity  on  the 
cavitation  bou  tdanes.  An  inevitable  question  is  which  of  the  abov  e 
two  deficiencit  s  s  more  severe? 

The  results  presented  in  Figs.  15  and  16  show  good  i  orrelat  ion  Ik* 
tween  the  "0**  <  avity  concept  and  the  results  of  Jakobson  and  Flo  berg 
fora  "universal  cavitation  parameter'*  of  -0.10.  Unfortunately,  thes< 
results  are  compared  with  "half  film"  solution  which  is  based  on  a 
cavitation  parameter  that  equals  zero.  A  more  realistic  comparison 
would  require  the  some  value  if  the  cavitation  parameter  for  all  th* 
solutions  compared  If  this  is  done,  the  solution  of  t  he*  cavitation  re 
gion  based  on  ‘quatiun  iSl  bec  omes 

6.|,,ia=f 

V  / 

where 

II  _ pmin  £  2 

Hence,  for  a  cnvitalion  parameter  -0,1  and  L/I)  ratios  0.25  and  0  5 
llmtn  becomes  - 1 .6  and  -0  4,  respectively.  It  is  clear  that  a  cavitation 
region  thus  obtained  from  (D-D  depends  on  /.//)  Therefore.  the 
eccentricity  Us  u*  and  load  capacity  are  also  L/l)  dependent  similar 
to  the  results  shown  in  Figs.  ]  » and  16  for  the  "0"  cavitat  ion. 

Additional  Reference 

17  Findlay,  b  A.,  "Ca*  Ration  in  Mechanical  Face  Seals."  Journal  ok 
f3tHRtCATK)Nl>:<  HNOLotiV.  V«  I.  90,  No  2.  Apr  1968 


N.  Tlpel4 

The  author  is  to  be  commended  for  this  paper  which  introduce*  m*w 
concepts  for  a  l»ettcr  underManding  of  cavitation  phenomena  m 
bearings.  The  initial  limitations,  regarding  the  short  Itearmg  ap 
proximation  and  the  submerged  arrangement  of  the  journal  Iwaring. 
only  partially  affect  the  validity  of  the  analysts  which  should  In  ex 
tended  to  other  conditions. 

The  previous  work  of  (ieoffrey  Taylor.  Pearson  and  Pitts  .ind 
(Ireiller  has  provided  up  in  a  certain  |toint,  valuable  information  on 
the  generation  and  on  the  number  of  streamers  in  sultcav  Rated  re 
gions.  However,  the  present  analysis  yields  specific  and  detailed  cal 
dilations,  to  lx*  applied  in  lubrication  problems. 

This  discusser  would  like  the  author's  comments  on  the  following 
points: 

1.  How  much  of  an  effect  does  the  surface  tension  have  on  the 
boundary  conditions  at  the  liquid  gas  interface? 

2.  The  usual  nonslip  boundary  conditions  let  a  thin  film  ol  In 
bricant  to  be  carried  by  the  journal,  even  in  sultcavilation  regions,  van 
the  assumption  of  dry  surfaces  be  critical  in  some  cases’’ 

3.  For  "I”  cavitation  the  diagrams.  Fig.  14,  do  not  show  value*  tor 
(  >  0.8.  In  opposition  to  short  fearing  theory,  should  the  ore  under 
ential  pressure  gradient  also  Ire  considered  at  higher  <  va  nes” 

In  conclusion,  the  preseni  work  is  a  significant  contribution  tu  a  still 
not  well  understood  chapter  of  lubrication.  The  extension  of  thi* 
analysis  to  other  types  of  liearings  should  Ire  welcome. 

Author’s  Closure 

The  author  is  grateful  to  Dr.  Ktsion  s  addition  ol  a  reference  to  t hi- 
earlier  work  of  Findlay,  who  was  concerned  with  the  lift  and  leakage 
of  a  lace  seal  as  affected  by  waviness  induced  cavitation,  it  nuiv  be 
mentioned  that  the  tracking  dynamic's  of  face  and  shall  seal*  would 
Ik*  dejrendent  on  the  cavitation  process. 

'The  alternate  solution  suggested  bv  Dr.  F.tsion  is  in  essence  an 
adaptation  ol  the <•  umbel  approximation  |2|.  which  is  at  best  a  con 
venient  short  out.  1'he  author  can  see  mi  merit  tn  of  fering  it  »s  a  li  xeb 
improvement.  1'he  crucial  concepts  to  lie  emphasized  hen  are  that 
1 1 )  1'he  Swift  -St  ieber  condition  |4,  5|  is  reconcilable  with  the  slioii 
I  waring  approximation  and  thus  should  he  asserted;  and  12)  Unle>* 
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I ht.*  ga,>  in  ti  led  with  lubricant.  no  pressure  ran  In*  developed.  As  lit 
the  ins ■wnl  pressure discontinuity at  th»-  Hint  re-establishment  point.* 
ti  litre  nnl\  hcs.ndtlinl  it  is  inevitable  to  allow  such  a  solution  if 'mu* 
.uufiN  I  In*  physical  realitv  n|  a  short  partial  an*  journal  lira  ring. 
<  'I* . trl . .  one  must  again  rek  on  matched  usvntpfnfir  expansion  |lt»| 
to  iriuoM*  hr  pressure  discont  miutv  at  both  entrance  and  exit 
Wilt'S. 

( *:u  ran  examine  i  hr  possible  merits  of  equation  il)- 11  as  a  partial 
mipro  cinent  over  the  eonvehlional  short  bearing  approximation. 
I'hiN  idea  hi  d  already  been  tried  in  a  study  «»t  squeezed  i  I  in  dampers 
I  Idj.  V  hde  •  here  was  a  definite  improvement  over  the  conventional 
-h  rt  bearing  approximation,  there  remained  significant  inexplicable 
dis  re  ianc  ies  from  ex|»erimental  lattN  One  inherent  fault  in  the 
(huith"l  Kl  inii  proposition  is  that  the  tavilation  domain  would  be 
uh\  i  V'  coot  lined  within  the  divergent  port  ion  ol  the  gap.  Figure  It) 
i  It  arlv  show  that  ihe  reestablishment  |»omi  inav  lie  in  the  eonvergent 
pt  int  it  tlu  gap.  depending  on  the  precise  value  ot  <  twhich  fixes 
: ,  For  i  istance.  it  the  numerical  value  ol  7rmtn  should  he  quite 
sm  II.  the  t  umhel  Ktsion  approach  should  yield  resiilis  very  dose 
ttit  ne  ull  I  dm  approximation,  which  are  contrndieted  bv  the  author’s 
? e -  tit'  h t  llM111(  =  — 0.01  Thus,  one  inav  salelv  conclude  that,  with 
wli.ie*  er  partial  improvement  may  he  ottered  by  the  use  of  equation 

•  I'  1  >.  some  uneerlaintv  associated  with  the  precise  location  of  the 
fill.,  n  establishment  point  will  result  The  argument  can  only  In* 
re». .l\e<l  by  a  subjective  cost  benefit  judgment. 

I  he  author  is  flattered  to  receive  the  kind  remarks  from  Dr.  Tipci. 
whose  contributions  to  the  field  ol  lubrication  theory  have  been  the 
noi  I've  ot  in  spiral  ion  of  l  he  author’s  ovv  n  el  forts.  Specific  quest  ions 
i.i  sed  bv  1  h  .  Tipci  .ire  answered  as  follows 

1  Warning  that  perfect  wetting  takes  place  between  the  luhri 
«.  »ii  *  a  id  i h.>  bearing  wall,  then  the  numscus  may  have  a  radius  of 
cnr.alure  which  is  lower  bounded  by  one  fourth  of  the  gap  where 

*  ,iv  ta  ini  breakup  takes  place.  <  hie  immediate  <  oiisequence  is  t hat 


flu*  film  pressure  at  tie  meniscus  ean  I**  I  »wt  r  !  Ian  the  cavity  pressure 
by  the  amount  of  |a//i,  where  a  is  the  surface  tenst<  n  of  the  lubri¬ 
cant  void  interlace.  Typically,  n  has  a  value  of  about  :  (i  dynes/cm. 
Thus,  il  t  he  gap  is  *J,a  1 1/ 10  mil),  then  the  associate*  1  pressure  jump 
al  I  be  meniscus  may  approach  tUJ  X  Id*'  dvnes/i  ji  -  ij.74  psi). 
However,  pressure  jump  is  only  one  lai  n|  el  the  tola  |  henonienoll. 
Al  the  breakup  meniscus,  using  the  )nc;'l  gap  as  the g-ornelric  scale, 
the  lubrication  theory  provides  the  up  t ream  lar-field  asymptotic 
limit  On  the  downstream  side,  the  flow  (assuming  absence  of 
si  lean  ers)  ultimately  becomes  a  thin  lilm  free  surlace  How  over  a 
moving  wall.  Transition  Ik* tween  these  two  vastly  di Here  it  flow  fields 
takes  plat  e  in  ihe  vicinity  of  the  menisc  is  [  Hi}:  itvrti.i »  fleets,  in  ad¬ 
dition  to  surface  tension,  may  lie  important  Iw  cause  di  *;ir  stress  el¬ 
icits  a  *e  minimal  on  the  downstream  skb  Instability  ot  lie  meniscus 
tli  w  ti  ‘Id  with  respect  to  a  transverse  w  a vt  is  |>r«  bal  Iv  the  cause  of 
the  si  reamer  structure,  Consequently.  he  pri  cise  pl.u  form  of  ihe 
meniscus  iMtundary  cannot  be  resolved  to  a  scale  of  t  lu  bearing  gap 
and  in  iv  further  be  altered  by  the  streamer  strut  lure  v.hrch  should 
have  a  wave  length  somewhere  between  the  bearing  gap  and  the 
bl  aring  lengl h  There  is  also  reason  to  believe  t  hat  the  y  wilt  St  ielier 
condition  should  tie  corrected  for  surface  tension  eflee1  *. 

J.  The  nonslip  lioundary  conditions  can  be  reconc  led  with  the 
“dry  cavity”  model  (associated  with  the  “I”  conlignr;  (ion I  by  ac- 
ce.ittng  wetness  to  In*  satisfied  by  a  film  of  molectil  ir  dimensions. 
Inerlic  effects,  however,  can  Ik*  responsible  for  residua  films  in  the 
“drv  cavil  v.” 

t  t 'ircunilerent ial  pressure  gradii  nt  is  propor  u  nnl  to  H~  \ 

Therefore,  given  a  small  but  finite  /./D,  the  short  healing  approxi¬ 
mation  would  eventually  become  inadoqu  ite  if  t  is  al  owed  to  ap- 
pr  ».hI  unity.  The  author  would  recoir  me  id  ant  ion  ;  gainst  using 
short  hearing  lesultsat  very  large  value  s  of  t  \  value  o'- 0.8  is  a  rea¬ 
sonable  upper  Ivound  for  trusting  short  hearing  results. 
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Cavitation  in  a  Short  Bearing  With 
Pressurized  Lubricant  Supply 

This  work  is  a  generalization  of  previous  treatments  of  the  cavitation  problem  for 
short  journal  bearings  [L  2\.  Cavitation  incipience  and  the  cavitation  pattern  are 
dependent  on  both  supply  and  cavitation  pressure  levels.  The  geometrical  pattern  of 
the  cavity  boundary  is  found  to  satisfy  a  universal  solution  which  depenas  uni  pi  civ 
on  the  cavitation  incipience  eccentricity  and  the  operating  eccentricity.  The  same 
eccentricity  parameters  also  uniquely  determine  the  friction  and  through-flow 
characteristics.  Radial  and  tangential  force  components  are  dependent  on  de 
supply /cavitation  pressure  ratio  in  addition  to  the  eccentricity  parameters  l  i  e 
results  are  applicable  to  journal  bearings,  shaft -seals,  and  squeeze-film  dampers. 


Introduction 

The  short  hearing  analysis  of  Dubois  and  Ocvirk  p]  is  a 
rare  combination  of  extreme  simplicity  and  near  validity.  I  or 
over  a  qua-  ter  of  a  centi  ry.  lubrication  engineers  have  used 
the  short  bearing  formulas  to  estimate  load  capacity  and 
friction  los*  of  plain  journal  bearings  It  has  been  as  common 
notion  that  e'er  for  not-so-short  beatings,  these  formulas  arc 
fair  approximations.  Although  numerical  computation  has 
become  standard  engineering  practice,  the  short  bearing 
analysis  remains  a  favoiite  approach  not  only  for  bearing 
designs  but  also  for  studying  rotor-bearing  dynamics,  and  for 
design  evaluations  of  shall  seals  and  dampers. 

It  is  a  common  practice  that  the  short  bearing  analysis  is 
performed  with  the  additional  approximation  of  retaining 
only  the  co  turned  film  pressure  in  the  converging  port. on  of 
the  bearing  gap  (ff-film).  This  turns  out  to  be  a  reasonable 
assumption  in  most  noi-so-short  bearing  designs  because 
lubricant  is  usually  copiously  supplied  with  the  feed  port 
located  in  the  bearing  wah  somewhere  in  the  divergent  portion 
of  the  bearing  gap.  The  ate  of  lubricant  flow  depends  on  the 
reservoir  pressure  level  and  the  flow  impedance  upstream  of 
the  feed  po  t.  The  reservoir  pressure  is  nominally  maintained 
above  the  jrrbient.  A  meniscus  boundary  is  formed  around 
the  feed  por:.  liven  though  the  miniscus  boundaty  .s  not 
precisely  the  maximum  gap  location,  the  overall  film  pressure 
distribution  docs  not  deviate  substantially  from  the  t  film 
approximation.  Ironicallv,  for  a  truly  short  configuration  (as 
may  be  in  the  case  of  a  journal  bearing,  a  shaft  seal,  or  a 
damper),  the  realises  concerning  lubricant  supply  and  the 
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film  extent  need  more  attention,  fypicallv ,  there  is  m»  teed 
port  in  the  bearing  wall.  The  lubricant  can  reach  the  yap  at 
either  end  and  can  enter  to  form  a  thin  film  only  it  the  e  :s  a 
region  where  the  film  pressure  i>  lower.  Consequent!'  ,  the 
lowest  pressure  level  which  can  be  sustained  by  the  luh-u.mt 
film  becomes  a  crucial  issue. 

Unfortunately  there  is  no  simple  and  practical  w.i\  ;<> 
predict  the  cavitation  incipience  pressure  in  fluid  Mm  d-'viccs 
despite  continued  interests  in  the  profession  (4,  5;.  The  best 
estimate  for  the  cavitation  incipience  pressure  is  some  a  here 
between  absolute  7ero  and  the  discharge  ambient.  In  the 
present  work,  the  cavitation  incipience  pressure  is  regarded  o 
be  a  given  constant,  which  not  only  decides  where  the  n  pure 
of  the  lubricant  film  to  take  place  but  also  gov  err  *  the 
location  of  the  full  film  formation  point  and  thus  the  ext  :nt  of 
cavitation. 

The  submerged  short  journal  bearing  represents  a  smirk* 
theoretical  model  to  allow  a  consistent  How  balance  bc  \wni 
lubricant  flow  into  the  divergent  gap  and  the  end  leakage  flow 
from  the  convergent  gap.  This  problem  was  thorough!) 
treated  in  a  receid  woik  (1).  It  was  found  that  if  a  uthmerced 
journal  bearing  is  to  carry  any  load,  the  cavitation  pressure 
must  be  suban!  knt  and  that  the  extent  of  cavitation  de  viuls 
on  both  the  cavitation  pressure  and  the  operating  ecccnti  icuv . 
Comparison  with  the  numerical  results  given  by  Jako^on 
and  Flobcrg  for  a  submerged  journal  bearing  of  finite  length 
(6)  demonstrated  that  the  improved  short  bearing  analysis  is 
valid  for  not-so-short  configurations.  According  to  the  same 
analysis,  the  submerged  short  bearing  has  been  found  to  be 
vulnerable  to  a  condition  of  lubricant  depletion  if  voids  it  the 
ambient  pressure  are  allowed  to  aggregate  near  ne  btai  tig 
ends. 

The  submerged  journal  bearing  configuration  ma>  not  be  a 
realistic  model.  Someomcs  lubricant  is  delivered  to  one  end  of 
the  bearing  at  a  more  or  less  elevated  pressure.  As  a  whole  the 
lubricant  passes  through  the  bearing  gap  and  is  discharged  to 
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the  ambient  conditon  at  the  other  end.  The  discharge  end  may 
or  may  not  be  flooded.  Shaft  seals  operate  precisely  in  this 
manner.  It  is  also  a  likely  confifuration  for  squeeze- film 
dampers.  The  elevated  feed  pressure  would  postpone 
cavitation  incipience  to  a  higher  eccentricity  and  would  reduce 
the  extent  of  cavitation  at  a  given  eccentricity  as  compared  to 
that  of  an  equivalent  submerged  bearing.  The  flow 
requirement  of  a  short  bearing  fed  from  one  end  was  studied 
by  Wakuri,  et  al.,  12);  an  ambient  cavity  was  assumed  to  exist 
at  the  discharge  end  of  the  short  bearing.  In  the  case  of  an 
internal  cavity,  the  feed  pressure  directly  influences  the  extent 
of  cavitation;  consequently,  the  through  flow  impedance 
varies  accordingly.  The  present  work  is  a  direct  extension  of 
earlier  studies  (I,  2)  by  allowing  the  pressure  at  the  feeding 
end  to  become  elevated,  by  including  the  possibility  of  a 
flooded  discharge  end,  and  by  rounding  out  the  results  to 
include  film  force  components  and  bearing  friction  torque 
together  with  the  through  flow  impedance. 

Analysis 

General  Short  Bearing  Theory.  Lubrication  theory  for  a 
short  journal  bearing  with  an  incompressible,  isoviscous  fluid 
obeys  the  following  different  ial  equation:  [3) 


=  H  3 


a// 

dr\ 


t  sin  17 

(1  -  (  cost*)1 


(3) 


(n0,n,)  are  invariant  with  respect  to  f  and  are  determined 
from  suitable  boundary  conditions. 

This  formulation  is  also  applicable  to  the  squeeze-f  Irr 
damper  undergoing  a  circular  whirl  motion.  The  coordinates 
are  fixed  to  the  whirl  motion  with  the  angular  coordinate 
aimed  against  the  sense  of  whirl.  u>  should  be  replaced  by  2r, 
twice  the  whirl  frequency,  and  *  now  represents  the  ratio  of 
the  orbit  radius  to  the  damper  clearance. 

Full  Width  Film  With  Circumferential  Feeding.  If  there  is 
no  cavitation  and  one  end  of  the  bearing  is  maintained  at 
Hs  >0,  then  upon  satisfying  the  boundary  conditions 

n(r=  -i)=ns;  n(r=n=o 

the  film  pressure  is  determined  by  a  special  form  of  equation 

(2V. 


n« Jo-mu-  (4) 

The  location  of  the  minimum  pressure  is  of  special  interest 
in  consideration  of  cavitation  since  the  value  of  the  minimum 
film  pressure  determines  whether  or  not  cavitation  would  take 
place.  The  minimum  pressure  is  located  by  the  conditions 


where 


6Mu >R2(L/D)2 

is  the  dimensionless  film  pressure  measured  above  ambient . 
For  an  axially  uniform  film,  H  »  I  -e  cost*,  equation  (1)  has 
the  general  solution 

f2 

n=n0  +  jHi  +  yCK»i>  (2) 

where 


an  1  ,  ao 

(5) 

dn  1  „ 

°~  “ar  2Ils+i"'<? 

(6) 

Equation  (5)  shows  that  the  circumferential  location  of  he 
minimum  pressure  is  not  dependent  on  ns  Equation  (61 
shows,  however,  that  unless  the  condition 

(7) 
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bO 
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i 
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ss 

radial  and  tangential 
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0 

= 
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force  components 
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/ 
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Z 
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cavity 
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w 
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€ 

r 

SI 

S3 

eccentricity  ratio,  e/C 
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Superscripts 

coordinate,  2 z/L 

nr 
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ac 

universal  variable 

tm 
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* 
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location  of  minimum 
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ns 
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vO 
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is  satisfied,  a  true  pressure  minimum  would  not  he  pre>ent  in 
the  film.  The  minimum  pressure  location  is  always  on  the 
discharge  side  in  the  axial  direction  =  n*/(2(?)  >  0).  1  he 
special  condition  Q  =  IIs/2  causes  the  minimum  pressure  to 
be  located  at  the  exit  edge;  i.e.t  fm  =  1 . 

Forms  of  Cavitation.  For  a  submerged  bearing,  previous 
work  indicated  the  existence  of  two  distinct  cavitation  pat¬ 
terns  [I].  In  the  “0"  cavitation  pattern,  the  pressure  in  the 
cavitated  legion  is  maintained  at  n( ;  film  breakup  ;ip- 
proaches  tie  cavitation  boundary  with  the  modified  Swift- 
Stieber  boandary  conditions  (7,  8).  In  the  *T*  cavitation 
pattern,  the  minimum  pressure  nc  exists  at  an  internal  point; 
cavitation  i egions  are  connected  to  the  ambient  edges  and  thus 
are  maintained  at  the  ambient  pressure.  They  are  in  effect 
“dry  pockets,** 

With  one  edge  of  the  bearing  maintained  at  an  elevated 
pressure,  there  are  similarly  these  two  major  varieties  of 
cavitation  orms.  Due  to  the  elevated  pressure  at  the  supply 
edge,  however,  symmetry  would  be  disrupted.  If  one  makes 
the  reasonable  assumption  that  only  the  liquid  lubricant  flows 
through  the  feed  line,  then  an  edge  cavity  can  only  be  present 
on  the  exit  (ambient)  edge.  The  entire  periphery  of  the 
pressurized  edge  is  thus  fully  effective  as  the  inlet  passage  for 
lubricant  t  irough  flow.  The  exit  cavity  may  be  either  wet, 
such  that  the  film  approaches  the  cavitation  boundary  with 
the  Swift-Stieber  condition;  or  it  may  be  dry,  and  the  ambient 
cavity  borders  a  film  with  some  subambient  region.  The 
solution  of  the  wet  exit  cavity  turns  out  to  be  a  special  case  of 
the  intermit  cavitation.  The  dry  exit  cavitation  problem 
requires  separate  attention. 

If  the  exit  cavity  is  dry,  the  film  pressure  near  the  cavity 
boundary  i:»  necessarily  subambient.  Therefore,  there  must  be 
at  least  a  lccal  region  along  the  exit  edge,  where  the  lubricant 
flows  back  from  the  “exit”  edge  into  the  film.  At  the  same 
time,  in  the  immediate  neighborhood,  a  dry  cavity  interrupts 
the  film.  Such  a  special  condition  is  possible  if  the  lubricant 
level  is  maintained  precisely  at  the  initiation  point  of  the  dry 
cavity.  Clearly,  this  is  an  unlikely  situation  except  possibly  in 
a  deliberate  laboratory  setup;  hence  no  further  attention  will 
be  given  to  the  dry  exit  cavity  problem. 

With  a  wet  exit  cavity,  while  a  zero  pressure  gradient 
condition  «s  satisfied  at  the  breakup  boundary,  the  film 
pressure  is  everywhere  higher  than  the  ambient.  Thus,  the 
lubricant  alw  ays  flows  out  along  that  portion  of  the  exit  edge 
which  is  reached  by  the  fluid  film.  At  the  same  time,  the  entire 
exit  edge  remains  at  the  ambient  pressure  including  that 
portion  which  borders  the  wet  cavity.  This  is  a  '‘flow¬ 
through”  condition  typical  of  most  bearing  and  seal  in¬ 
stallations;  the  film  pressure  is  everywhere  higher  than  (he 
ambient  even  if  the  lubricant  is  physically  capable  of  a 
subambient  >tate.  The  exit  edge  can  be  flooded  by  the 
lubricant  (lor  example  by  restricting  the  exit  flow),  then  the 
cavity  will  be  located  inboard  and  the  pressure  in  the  cavity 
should  be  approximately  it  the  effective  vapor  pressure  of  the 
lubricant  and  can  be  substantially  subambient. 


Cavitaticn  Incipience.  The  full  film  solution,  equation  (4), 
is  applicable  so  long  as  the  minimum  film  pressure  is  higher 
than  the  cavitation  pressure;  i,e.f  II  >  Ilr.  The  condition  of 
cavitation  incipience  is  reached  when  the  minimum  film 
pressure  is  exactly  equal  to  the  cavitation  pressure. 

According  to  equation  (5),  the  circumferential  location  of 
the  minimum  piessure  point  (in  a  full  film)  is  given  by  *tQ/drj 
*  0.  Since  this  condition  does  not  depend  on  the  values  of  \\s 
and  nr,  the  relationship  between  the  circumferential  location 
of  cavitation  incipience,  tj(,  and  the  incipience  eccentricity 
ratio,  as  determined  in  the  previous  study  (2)  remains 
valid.  The  »alue  of  r,t  however,  is  now  dependent  on  both  IIV 


Fig.  1  Eccentricity  ratio  tor  cavitation  incipience 


J.“  <1.1  .  >.  j  .<  .  v 1  It.  r  *  '  1 


Fig.  2  Supply  pressure  level 


and  nr.  Upon  setting  Q  =  Q,  =  Qfi,,  r;,)  and  s',,.  =  $*  in 
equation  (6),  one  obtains 

Us  =  2Q,  f,  (8) 

Then  substituting  into  equation  (4)  aJong  with  (II,  ",  Q ), 
respectively,  identified  as  (Il(  ,  f, ,  Q{ ),  it  is  found 

-ii<  =<?id  -r,):/2  (*) 

Since  ij,  depends  uniquely  on  e,,  Q ,  can  be  regarded  as  a 
function  of  t,  also.  Accordingly,  equations  (8)  and  (9)  can  be 
regarded  as  nA  (*,,  £)and  -  nr  («,,  *",). 

Physically,  it  is  more  natural  to  regard  (*,,  f,)  as  functions 
of  (-Hc,  II*).  To  obtain  such  relationships,  one  car  first 
eliminate  Q,  between  equations  (8)  and  (9),  then,  after  some 
rearrangements,  it  is  found 

f1  =  i+2(-n1 7iis)[i*vT4(-ri$/nc)i  do) 

Then,  substituting  into  equation  (8),  one  obtains 

<?,  =<nv/2)M-n()W(-n<)[(-n(  >  +  n,]  on 

Since  Q,  is  a  unique  function  of  the  inverse  would  be 
true,  so  that  equation  (II)  yields  indirectly  the  function  ot  c, 
( -II< ,  Ilv).  Actual  computation  of «,((?,).  however,  requires 
a  numerical  procedure.  Graphical  representation  of  i,  (  -  II,  , 
nA)  is  given  in  Fig.  1. 

<-nr,  ns)  and  <£,  t ()  are  alternative  parametric  sets  to 
specify  the  condition  of  cavitation  incipience  Fig.  2  together 
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wiih  Fig.  1  furnish  the  complete  relationships  among  these 
parameters  in  the  graphical  form.  In  subsequent  analysis  to 
treat  the  captation  problem,  the  condition  t  >  (,  is  of  in¬ 
terest.  Thus  the  parametric  set  (£,<,)  will  be  assumed  to  be 
given. 

Universal  Solution  of  the  Cavitation  Region.  For  any 

combination  cf  a  Finite  and  positive  (nonvanishing)  value  of 
-  Il<  together  with  a  finite  and  non-negative  value  of  ns,  the 
axial  locatioi  of  the  incipience  point  is  internal;  i.e.,  0  <  £  < 
1.0. 

The  case  (Tv  =  0,  f,  =0)  corresponds  to  the  “0”  cavitation 
prolem  of  a  submerged  journal  which  was  treated  previously 
(1).  If  the  eccentricity  ratio  exceeds  the  incipience  value,  t  > 
t,  ,  then  the  cavitation  pressure  nr  would  be  reached  at  0  < 
»h)(  <  V  Th-  subscript  “W*  designates  the  condition  of 
“breakup  in  nation.”  Since  the  breakup  initiation  should  also 
be  at  the  axial  minimum  of  the  film  pressure,  equation  (6) 
must  apply.  Substituting  (II<  ,  fhl,  Qh, ),  respectively  for  (II,  F, 
Q)  in  equation  (4),  and  (jy  ,  0^),  respectively,  for  (£*,„,  Q)  in 
equation  (6).  one  obtains 

n,  ■=  -Jd  (121 

Us  =2{btQbt  (1^) 

Eliminatirg  -  II s ,  one  obtains 

-iv  ■•(?*,(  i-r^)2/2  (U) 

Comparing  equations  (13)  and  (14)  with  equations  (8)  and 

(9),  one  concludes  that  since  (fA/>  Qh)  have  the  same  func¬ 
tional  relations  of  ( -  IIr,  Ils)  as  (f,,  Qt)f  it  follows 

f*-  ft;  Q»=Q,  UM 

Noting  that  f,  depends  only  on  the  parameter  (-Il(  /lls) 
according  to  equation  (10>,  the  value  of  Jy,  =  f(  is  a  direct 
indicator  of  me  supply  pressure  (relative  to  the  cavitation 
pressure).  rjh  ,  on  the  other  hand,  can  be  solved  from  equation 
(3)  upon  substituting  Qh  =  Q  (f,)  for  Q ,  and  therefore  is 
dependent  only  on  (<, «,)  in  the  explicit  sense.  Computation  of 
!)„,((,  <,)•  however,  has  to  be  performed  numerically  since  the 
governing  re  aiionship  is  transcendental. 

t  si  irj^ -<  cosi^,)’  =0  (16) 

The  break  up  process  continues  beyond  gA,.  An  internal 
cavitation  domain  separates  the  two  “side  bands”  which 
border  either  end  of  the  bearing.  The  pressure  field  in  the 
sideband  is  i  mvcrsal  because  the  pressure  of  either  end  of  the 
bearing  is  cjr slant,  whereas  the  Swift-Stieber  condition  is 
applicable  a:  he  cavitation  border.  The  universal  pressure 
profile  can  thus  be  written  as 

iWe./o  a?) 

(II,  f)  are,  tespcctivcly,  the  universal  film  pressure  and  the 
universal  axial  coordinate  for  either  sideband.  They  arc 
linearly  related  to  (11.  f)  with  appropriate  scale  factors  and 
offsets;  i.e., 

riMH-IU/AJl;  f=(jv-r0>/^f  (18) 

The  scale  factors  and  the  offsets  are  listed  in  Table  1 . 

In  the  uni’  ersal  representation,  as  shown  by  equation  (17), 
the  width  ot  the  sideband  depends  parametrically  on  (< .  O 
and  is  a  function  of  i?  It  assumes  its  minimum  value  when  Q 
reaches  its  ma*  imum  at 

9.  -  cos  1 1  (Vl  f  24«J  -  l)/(4<)|  (19) 

The  subscript  “w”  refers  to  the  “waisl“  of  the  sidehard.  If 
cavitation  hid  not  taken  place,  would  locate  the  minimum 
film  pressure;  i)„  is  functionally  dependent  on  c  identically  as 
i),  is  on  i,. 


Fig.  3  Geometrical  perimeters  ot  the  universal  cavitation  pattern 


Beyond  n*.  the  sideband  widens  with  rj  as  the  axial  How, 
drawn  in  from  the  end  adds  to  the  circumferentially  tians- 
ported  flow,  through  the  internal  cavity,  to  cause  the  film  to 
fill  back.  The  universal  pressure  profile  is  given  by 

II*  “(f-ff)/fr  +  «f-frKO/e,)  (20) 

{(  is  the  cavitation  boundary  in  the  “fill-back”  region.  Note 
that  the  boundary  values  of  ri  at  $  -  (0,  f,)  are  satisfied, 
whereas  the  Swift-Stieber  condition  is  no  longer  enforced 
except  at  rj„  where  f,  =  jcH . 

At  the  fill-back  boundary,  flow  transported  through  the 
cavity  from  ij(ft,  =  jp, ) .  which  is  given  by  equation  (17), 
combines  with  the  axial  in-flow  through  the  widening 
sideband  to  fill  the  gap  of  the  reformed  film.  For  the  *hort 
bearing  approximation,  only  the  Couette  component  is 
retained  in  the  circumferential  flow  of  the  reformed  film.  The 
relation  between  f,  and  >  tj*  is  same  as  that  for  the  sub¬ 
merged  bearing.  Upon  approximating  the  local,  transported 
flow  by  the  average  value  (l),  and  adopting  the  universal 
representation  the  formula  of  interest  is 

( H - Hw)jyJ  =(».-«*,)<?,/ 0.  +  in i»)  -  F(  >  10,  (2D 
where 


AM- Jm** 


sin 

12 


=  ('  4  j  ,  ,'(l  *  4  )sin''  4  4  ,  si,'-7 


(22) 

If 

cosi|w  =  - — ;  cosij  (2-1 1 

(If-)  J  -V 

and 

Hm  =  I  -t  cos**  (24| 

The  Fill-back  process  is  terminated  when  Jy  becomes  unity. 
The  universal  form  of  jy  and  jy  ensures  that  the  internal 
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Fig.  4  Cavitation  patterns 


cavitation  begms  and  ends  together.  For  convenience  of 
future  reference,  the  tetminaiion  location  is  designated  tj,  . 
The  computations  previously  performed  for  the  "0” 
cavitation  problem  of  the  submerged  bearing  (1]  are  in  fact 
directly  applicable  to  the  internal  cavitation  problem  viewed 
in  the  universal  form. 

In  the  universal  representation,  geometrical  parameters  of 
the  internal  cavitation  pattern  can  be  summarized  as  the 
circumferential  locations  of  the  breakup  initiation  point,  the 
film  waist,  and  the  cavitation  termination  point  together  with 
the  width  fraction  of  the  film  waist.  They  are  uniquely  cefined 
by  («,, c)  as  graphically  shown  in  Fig.  3. 


Table  1  Scale  factors  and  offsets  in  side  bands 


Side  band 

Inlet 

Exit 

Pressure  oft  scr  n0 

»< 

n<~ 

Pressure  scale  All 

Us  IV 

-ii, 

Axial  offset  f0 

-  1 

1 

Axial  scale  Af 

1,-1 

f,  i 

Results 

Cavltallnn  Patterns.  The  pressure  parameters  (  ~  Tl(  ,  II  v) 


affect  the  cavitation  pattern  in  two  ways.  F irstly  thev 
determine  t,  according  to  Fig.  I.  It  may  be  noted  that  the 
lower  bound  of  * ,  depends  only  on  II s  as  given  by  the  line  f(  = 
1.0  in  Fig.  2.  <,  together  with  the  operating  eccentricity  <  fix 
the  circumferential  characteristics  of  the  cavitation  boundary; 
i.e.,  the  locations  of  breakup  initiation,  the  waist,  a  id  the 
termination  point  along  with  the  waist  width  ratio  £„  as 
shown  in  Fig.  3.  Secondly,  the  two  "halves”  of  the  cavitation 
pattern  divide  at  £  being  a  function  of  the  ratio  ( -  llr/Hs ) 
as  given  by  equation  (10),  and  (1  +  £,  I  -  £)  are,  tespec 
lively,  width  scales  of  the  two  halves  in  the  universal  sol  Jtion. 

These  effects  are  shown  in  Fig.  4.  Figure  4  (<r)  corresponds 
to  the  condition  («,  =  0.025,  r  =  0.10).  Foui  cavitation 
patterns  are  shown.  Each  plan  form  covers  the  entire 
developed  film;  the  lower  and  upper  boundaries  are  the 
minimum  gap  locations.  From  left  to  right  the  captation 
patterns  correspond  in  sequence  the  ratio  -Il(  ,T1S  =  (». 
1.0,  0.1,  0.0).  In  this  case,  because  f,  is  quite  small 
(representing  a  lubricant  with  an  effective  vapor  pressure  .er\ 
close  to  the  ambient  (and  a  very  low  feed  pressutc)  the 
cavitation  region  spans  about  2/3  of  the  total  periphery  wen 
(hough  the  operating  eccentricity  is  quite  modes'.  In  Fig. 
0),  with  f,  staying  at  0.025,  *  is  increased  to  0.3.  The  lour 
patterns  correspond  to  the  same  sequence  of  the  pressure 
ratio.  The  larger  operating  eccentricity  is  seen  to  be  ac 
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companies  by  the  enlargement  of  the  cavity  in  both  axial  and 
circumferenctial  directions.  Figure  4  (c)  corresponds  to  (*,  = 
0.2,  t  =  0.3).  The  m  >derately  large  t ,  requires  either  a 
relatively  low  vapor  pressure  lubricant  (about  10  times  more 
remote  from  the  ambient  than  the  previous  cases)  or  a 
moderately  high  feed  pressure  or  some  suitable  combinai  ton 
of  the  two  pressure  parameters.  Figure  4  (d)  shows  the 
cavitation  paterns  for  a  rather  large  operating  eccentricity. 
Note  that  the  length  of  the  cavity  remain  at  about  1/2  the  total 
periphery  while  its  width  has  increased  substantially. 

It  is  quite  clear  that  the  feed  pressure  is  an  effective  means 
to  limit  the  circumferential  extent  of  cavitation. 

Formulas  for  Global  Characteristics.  The  essential  per¬ 
formance  parameters  of  a  journal  bearing  include  the  force 
components, 


{WR,WT) 


(4M4) (i) 

*=  J  (cos?),  —  siiiif(di)  ^  Ud^ 
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the  friction  torque, 
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The  terms  associated  with  the  subscript 
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Because  the  integrands  take  on  different  functional  forms 
in  the  full-film  and  cavitation  regions,  it  is  convenient  to 
separate  the  integrals  into  parts  pertaining  to  the  two  regions. 
One  can  therefore  write 
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TC 
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*F  **  are  con¬ 
tributed  by  the  full-film  region,  which  spans  i),  -  2k s  17  ^ 
11^.  The  subscript  “C”  identifies  the  contribution  to  the 
cavitation  domain,  r)b,  s  i)  s  ij,  .  It  is  not  necessary  to 
separate  the  supply  flow'  into  two  pans  because  the  cavity 
blocks  axial  rhrough  flow  so  that  £  is  entirely  due  to  the  full- 
film  region. 

For  short  bearings,  contribution  to  friction  torque  by 
circumferential  pressure  gradient  is  of  O{(/70)2),  and  is 
commonly  neglected.  In  the  cavitated  region,  i)b,  s  i)  s  i)n 
the  friction  torque  (7)r  is  made  up  from  two  parts;  one  part  is 
due  to  the  continuous  films  between  either  end  of  the  bearing^ 
and  the  cavity  border  (0  <  f  <  {b  for  s  ij  <  ij*  and  0  <  f 
<  fy  for  ij,„  s  7 1  £  ip),  the  second  part  is  due  to  the  striated 
streamers  wirhin  the  cavity  (9).  For  computing  the  streamer 
contribution  to  (t)c •,  the  local  gap  at  film  breakup  is  used  in 
the  range  ij*  3  s  and  the  average  value  is  used  in 
the  range  if  ,  $))$>),, 

In  summary,  the  global  characteristics  in  the  dimensionless 
form  are 
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Various  subscripts  arc  used  to  designate  the  appropriate 
integration  limit.  One  may  note  that  T  and  i  can  be  fully 
described  in  erms  of  (<,,  r),  which  uniquely  determine  (i y„f,  ry, > 
as  well  as  jyij,  ,  s  ry  <  ?jw)and  fy(iyv  s  y  <,  ty,),  without 
further  specification  of  (II  v ,  flr).  ff'r),  however,  are 

influenced  by  (IIS,  n(  ),  in  addition  to  (c,  <4)  because  of  the 
pressurization  effects  both  in  the  full-film  region,  equations 
(29a)  and  (30er),  and  in  the  cavitation  region,  equations  (29c) 
and  (30c). 

Representative  Results.  Graphs  for  the  dimensionless 
friction  torqje  are  shown  in  Fig.  5.  A  value  of  unity  indicates 
agreement  with  Petrov  s  law  [10).  Upper  and  lower  bounds 
are  respective  full  film  and  e-film  frictions  with  allowance  for 
eccentricity.  It  is  observed  that  even  though  the  film  may  be 
extensively  capitated;  for  instance,  for  («,  =  0.025,  €  =  0.3), 
as  shown  by  the  cavitation  pattern  of  Fig.  4  (b),  the  torque  is 
substantially  higher  than  e  film  result.  The  effect  of  enlarging 
the  cavity,  however,  is  evident  for  ti  =  0.025,  as  the  curve 
dips  down  and  reaches  a  minimum  value  near  t  =  0.4. 
Competing  effects  of  eccentricity  and  cavitation  nullify  each 
other  such  that  near  agreement  with  Petrov's  taw  is  indicated 
for  (0.025  Sf,s  0.20,0  «;  «  s  0.55). 

Graphs  for  the  flow  function  are  shown  in  Fig.  6.  The  same 
results  were  previously  obtained  by  Wakuri,  et  ah,  [2]  with  the 
interpretation  restricted  to  the  case  of  an  ambient  exit  cavity, 
so  that  each  value  of  #,  is  associated  with  a  unique  value  of 
n5.  It  is  dear,  nevertheless,  that  these  results  remain  ap¬ 
plicable  to  the  bearing  with  a  flooded  exit  end,  such  that  an 
internal  cavity  would  sustain  a  subambient  pressure.  In  the 
latter  situation,  if  both  ^  and  <  are  available  from  test  data 
and  n,  is  known;  then  Fig.  6  allows  one  to  determine  t,,  and 
subsequently  Fig.  I  can  be  used  to  determine  -  Ilr  for  the 
particular  combination  of  (fl5,  c,).  In  this  manner,  a  j7ow- 
ecctn  tricity  experiment  would  become  a  basis  to  determine  the 
cavitation  pressure. 
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Fig.  T  Hydrostatic  offsets  on  radial  fores 
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R(t.  6  Hydrostatic  elfacta  on  tangantial  fore# 


roa:£5  AT  small  ICCENTaiClTY 
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Fig.9  Hydrostatic  «f  facts  on  tores  componsots  at  small  sccsntncity 
for*,  *  .02$ 

It  was  noted  previously  that  the  force  components  cannot 
be  uniquely  specified  in  terms  of  (i,  <  ).  This  is  because  the 
pressure  profile  in  the  cavitation  domain  takes  on  a  drastically 
different  character  than  thar  in  a  full-width  film.  In  a  full-film 
condition  as  shown  by  equation  (4),  the  contribution  of  I1N  is 
routionally  symmetrical  and  thus  would  not  result  in  a  net 
force.  The  lull-film  journal  bearing  force  is  entirely  due  to  the 
hydrodynamic  wedge  effect  represented  by  Q.  In  a  cavitaied 
region,  even  if  one  continues  to  measure  the  pressure  profile 


from  the  straight  line  which  connects  the  pressure  levels  of  the 
two  ends,  its  characterization  requires  representation  of  the 
pressure  and  width  scales  as  indicated  in  Table  1 .  The  roles  of 
(II* ,  - 11( )  in  the  pressure  scales  are  self-evident.  The  width 
scales  contain  the  parameter  f(,  which  in  turn  is  a  function  of 
the  pressure  ratio  ( -Flf  /IIS).  Since  e,  already  provides  one 
constraint  between  (IIS,  -hr),  a  second  parameter  is  suf¬ 
ficient  to  specify  both  (IIS,  -Ilr).  For  the  latter  purpose, 
three  values  of  the  pressure  ratio,  -n(  /nv  *  (a#,  0. 1 , 0),  arc 
used  to  examine  the  effects  on  the  forces  which  are  not  in¬ 
cluded  in  the  hydrodynamic  components  as  represented  by 
equations  (29b)  and  ( 30£).  The  same  conditions  were  usee*  to 
illustrate  the  sequence  of  successively  increasing  asymmetry  in 
the  cavitation  pattern  in  Figs.  4(a)  through  4(d).  In  Figs.  7  and 
8,  radial  and  tangential  forces  are,  respectively,  plotted 
against  <  in  two  groups,  each  of  which  includes  the  three 
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Fig.  1 1(a)  Dinxfiftiontos*  langintlil  force  •  submerged  bearing 


values  of  tie  pressure  ratio  for  a  given  The  two  values  of  t , 
are  (0.1, 0  4).  For  comparison,  the  hydrodynamic  component 
is  drawn  in  dashes.  Both  radial  and  tangential  forces  show 
moderate  spreads  among  various  values  of  the  pressure  ratio 
at  moderate  values  of  t  >  In  the  same  range  of  <,  the 
hydrodynamic  components  are  not  good  approximations  of 
the  actual  orces.  Further  scrutiny  of  this  question  at  small  («, , 
€)  is  referr  :d  to  Fig.  9.  So  far  as  the  radial  force  is  concerned 
(for  c  a  <  ),  the  hydrodynamic  component  is  lower  than  the 
full  force  at  a  small  t,  and  can  be  an  order  of  magnitude  coo 
large  at  moderately  large  values  of  e(.  The  hydrodynamic 
component  of  the  tangential  force  is  always  smaller  than  the 
full  value,  it  breaks  awa>  from  the  full-film  value  too  rapidly. 
The  two  conditions  (fl.s  =  0,  IIr  =  0)  represent  the  extreme 
conditions  of  symmetry  For  any  combination  of  fa,,  <),  all 
solutions  of  the  forces  are  somewhere  between  those 
corresponding  to  the  extreme  conditions.  Figures  10  (a)  and 
( b )  contain  curves  for  the  radial  farce  for  the  two  respective 
conditions  (IIS  =  0,  Il(  =  0).  Curves  for  the  tangential  force 
are  given  in  Figs.  1 1  (a)  and  ( b ).  It  turns  out  that  for  any 
finite,  combinations  of  (II5,  II( ),  the  following  empirical 
formula  yields  very  accurate  results: 

24.721  - ilcK »%, .(;,,VB0^ns(n#.,.Ilr  « 

(WV:n,.iif  ns  +  24.72(-n<  ) 

where  ft  nay  be  either  WR  or  ft'r. 

Conclusions 

Pressurization  of  lubricant  supply  at  one  end  of  a  short 
bearing  would  shift  the  cavitation  pattern  toward  the 
discharge  side.  For  the  same  cavitation  (subambient )  pressure, 
the  eccentricity  ratio  for  cavitation  incipience  would  be  in¬ 
creased  by  supply  pressurization. 


< 


Fig.  1 1(b)  Dlmanstonlaat  tangential  tore#  amblanl  ax  it  captation 

Size  parameters  of  the  cavitation  pattern  are  dependent  or 
fa,,  f)  only,  so  long  as  the  short  bearing  theory  is  valid  Axial 
location  of  the  cavity  and  the  degree  of  geometrical  asym 
metry  are  dependent  on  the  supply/cavitation  pressure  ratio 
(both  pressures  are  measured  above  the  ambient).  For  anv 
finite  «,  >  0,  there  is  an  upper  bound  of  the  extern  of 
cavitation  as  <  —  1.0. 

The  cavity  can  occupy  an  internal  domain  only  when  the 
exit  edge  is  flooded  by  lubricant  (with  negligible 
pressurization),  and  the  cavity  pressure  would  then  hi 
necessarily  subambient.  If  the  cavitation  pressure  is  preciseK 
equal  to  the  ambient  pressure  or  if  the  exit  lubricant  is  allowed 
to  drain  away  freely,  then  the  exit  edge  would  be  only  part  iall> 
connected  to  the  fluid  film  and  the  cavity  must  communicate 
to  the  ambient  directly.  The  drain  away  condition  enfoxes  ar 
ambient  condition  of  the  cavity  even  if  the  lubricant  can 
inherently  sustain  a  subambient  condition. 

Dimensionless  friction  and  through-flow  function  ;an  hi 
fully  defined  by  fa,,  f).  Competitive  effects  of  eccentric  tv  and 
cavitation  on  friction  results  in  the  near  agreement  wit! 
Pehrov’s  law  for  a  wide  range  of  conditions.  The  through 
flow  function  is  strongly  dependent  on  the  value  of  f  ir 
practical  operating  condition  of  f.  Thus,  validity  of  :hi 
present  analysis  can  be  verified  by  correlating  the  mutual 
dependence  between  the  cavitation  pressure  and  the  through 
flow  function  at  fixed  supply  pressure  and  the  operating 
eccentricity. 

Radial  and  tangential  force  components  (as  functions  of 
operating  eccentricity  <)  are  dependent  on  the  sup 
ply /cavitation  pressure  ratio  in  addition  to  the  determination 
of  the  cavitation  incipience  eccentricity  t,  by  the  supply  and 
cavitation  pressures.  An  accurate  empirical  formula  for  both 
force  components  has  been  found  to  permit  simple  in 
terpolation  from  their  respective  values  at  the  extreme  con 
ditions  of  zero  supply  pressure  and  zero  cavitation  pressure 
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DISCUSSION 


I.  Etsion2 

The  authors  have  presented  a  valuable  solution  to  the 
problem  of  short  journal  bearings  fed  from  either  one  or  two 
ends.  This  problem  is  also  of  great  interest  in  squeeze  film 
dampers  and  mechanical  tace  seals.  The  exact  solution  for  the 
cavitation  snape  and  bearing  performance,  as  derived  by  ihe 
authors,  is  quite  complex.  However,  such  exact  solutions  are 
essential  as  a  basis  for  assessment  of  less  accurate,  simpler 
solutions  that  may  be  used  for  practical  purposes. 

A  possible  simpler  solution,  for  example,  is  the  one  based 
on  equation  (4)  with  an  additional  condition  replacing  any  II 
<  n(  by  II  a*  nt  (the  half  Sommerfeld  condition).  The 
cavitation  shape  {*(17)  in  this  case  is  found  from  a  solution  of 
equation  (4  with  n  =  II, .  The  accurate  shape  of  the  cavity 
and,  more  important,  the  flow  continuity  on  the  cavitation 
boundaries  are  not  preserved  by  this  simple  solution. 
Nevertheless,  the  overall  results,  especially  for  friction  torque 
and  load  carrying  capacity,  may  be  fairly  close  to  these  ob¬ 
tained  from  the  much  more  complex,  exact  solution  presented 
in  this  paper.  The  authors  have  apparently  attempied  a 
comparison  between  their  exact  solution  and  available 
aproximations.  In  Fig.  5,  for  example,  the  results  for  the 
friction  torque  are  compared  with  these  obtained  for  eccentric 
full  film  and  eccentric  t  film.  The  latter  is  basically  the  half 
Sommerfeld  solution  for  the  case  -IIt  =n,  =0.  Un¬ 
fortunately,  the  r  film  results,  shown  in  Fig.  5,  do  not  ac¬ 
count  for  th:  striation  affect  on  the  friction  torque,  leading  to 
the  erronecus  conclusion  that  .  .  .the  torque  is  sub¬ 
stantially  higher  than  it  film  result.  .  .  Table  2  presents 
results  of  the  dimensionless  friction  torque  for  the  ir  film  with 
the  striation  effect  included.  These  results  were  obtained  by 
dividing  the  dimensionless  friction  force  /'of  reference  ( 1 1 1  by 
2t. 

Fm  the  results  in  Table  2  and  those  in  Fig.  5,  it  is  clear  that 
the  t  film,  or  rhe  half  Sommerfeld  solution  for  -  II,  =  n,  =  0, 
approximates  very  well  the  exact  solution  in  cases  of  suitable 
caviation  incipient  (e.g.  «,  «  0.025  and  e,  »  0.05).  The  exact 
solution  for  higher  values  of  t,  can  be  well  approximated  if 


more  realistic  (less  than  ambient)  cavitation  pressures  art  used 
in  the  half  Sommerfeld  solution.  Hence,  it  can  be  concluded 
that  the  half  Sommerfeld  solution  based  on  equation  (4) 
approximates  quite  well  the  exact  friction  torque. 

With  regard  to  the  lubricant  through  flow  it  can  be  snown 
that  the  half  Sommerfeld  solution  yields  the  same  results  as 
the  exact  solution.  This  is  due  to  the  fact  that  the  net 
hydrodynamic  component  of  the  flow  in  a  submerged  bearing 
is  zero  and  tjie  through  flow  is  exclusively  a  result  of  the 
hydrostatic  pressure  differential  across  the  bearing. 

It  would  be  interesting  to  compare  the  results  of  the  radial 
and  langenial  forces  obtained  from  (he  half  Sommerfeld 
solution  with  these  obtained  from  the  present  exact  sole  lion. 
Such  a  comparison  would  settle  the  yet  unresolved  issue  of  the 
accuracy  of  using  the  half  Sommerfeld  condition  for  ap¬ 
proximating  cavitation  in  narrow  configurations  like  face 
seals  and  squeeze  film  dampers. 

Another  problem,  which  was  not  given  much  attention  in 
the  past,  is  the  assumption  of  constant  pressure  inside  the 
cavitation  zone  and  on  its  boundaries.  This  assumption  which 
is  used  here  by  the  authors  and  was  used  extensively  by  other 
researchers  in  the  past  may  not  be  correct  for  enclosed  cavities 
such  as  in  submerged  journal  bearings.  In  a  recent  work  [12], 
where  the  pressure  field  in  the  cavitation  zone  of  a  submerged 
journal  bearing  was  maped,  it  was  found  that  significant 
pressure  variations  occur  inside  the  cavity.  In  cases  where  the 
full  film  formation  point  penetrates  into  the  converging 
portion  of  the  bearing  clearance,  high  cavity  pressures  were 
recorded  that  exceeded  even  the  elevated  pressure  at  the  ends 
of  the  bearing.  This  phenomenon  was  attributed  to  air 
solubility  of  the  lubricating  oil  which  probably  plays  an 
important  role  in  the  cavitation  mechanism. 
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Tabic  2  DlmeaaloalcM  friction  torque  for  a  film  with  striation  effect  included 

r  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8 

TCI^MR'ODL)  0.96  0.94  0.93  0.93  0.96  1.02  Ml _ 129 


0.9 

1.75 
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The  authors  concur  with  Dr.  Etsion  that  if  one  allows  for 
the  contribution  of  shear  from  striations  in  the  divergent 
portion  of  the  gap*  the  »-film  analysis  yields  rather 
satisfactory  results.  However,  one  may  justifiably  argue  that 
the  striation  correction  itself  implies  the  need  to  obey  overall 
continuity  of  flow;  otherwise,  as  indicated  by  Fig.  5,  the 
friction  torque  would  be  consistently  and  substantially  un- 
derestima  ed.  The  important  message  revealed  by  the  more 
cautious  calculation  is  that  the  simultaneous  effects  of  ec¬ 
centricity,  cavitation,  and  striation  tend  to  balance  one 
another  such  that  Petrov’s  law  is  quite  satisfactory  in  its 
simplest  form. 

It  is  quite  natural  to  query  to  what  extent  film  force 
calculation  may  be  adequately  approximated  by  the  ir-film 
results.  Th  s  question  was  examined  for  the  submerged 
journal  rearing  |1]  and  graphical  comparisons  were 
separately  presented  for  the  tangential  and  radial  force 
componer  ts.  These  graphs  are  reproduced  as  Figs.  A- 1  and  A- 
2.  If  one  modifies  the  Ocvirk-Gumbel  approximation  by 
requiring  fl  >  n,  with  a  suitably  selected  value  of  IT  <  0, 
then  the  corresponding  results  would  be  between  those  of  the 
standard  r-rilin  results  in,  =  0)  and  the  full  film  results  (Ilr 
<  nmin).  It  is  clear  that  since  some  of  the  curves  in  these 
Figures  arc  not  contained  between  these  “extreme  cases,”  the 
modified  Ocvirk-Gumbel  approximation  has  limited 
usefulness .  The  basic  fault  of  the  Gumbel  approximation  is  its 
inability  to  describe  the  boundary  of  the  re-established  Film, 
which  car  icach  into  the  convergent  portion  of  the  bearing 
gap.  This  conclusion  remains  valid  even  if  one  edge  of  the 
bearing  is  pressurized. 

Dr.  Etsion's  comments  regarding  the  hydrostatic  through 
flow  prompted  the  author  to  re-examine  the  role  of  cavitation 
in  this  respect.  Equation  (32)  and  Fig.  6  are  based  on  the  axial 
flow  in  the  full-film  (17,  -2»  <  rj  <  r)h<)  at  the  mid  plane 
=  0).  Wi  h  pressurizaiion  of  one  edge  of  the  bearing, 
however,  the  “cavitation  axis”  is  shifted  to  f,  >  0.  Ac¬ 
cordingly,  the  correct  hydrostatic  through  How  should  be 
calculated  as  the  axial  flow  in  the  full  film  at  £  and  equation 
(32)  should  be  corrected  by  a  “cavitation  leakage”  term: 
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Fig.  A 2  Radial-tore*  for  **0"  cavitation 


Fig.  A-3  Total  flow  function  equations  (32)  and  (A- 1 ) 


With  the  addition  of  “cavitation  leakage,”  the  total  fov* 
function  is  as  shown  in  Fig.  A-3. 

The  formal  simplicity  of  Equation  (A-l)  motivated  a  search 
for  similarly  simple  formulae  for  the  effects  of  pressurization 
on  the  force  components.  This  was  found  possible  upo 
recognizing 

IV -211  r«<l  +  f?K?(  (A-2i 

rn-d+JifK?,  <A-3> 

F  as  a  function  of  (IV,  -  H( )  can  be  readily  calculated  with 
the  aid  of  Equation  (10).  Therefore,  in  place  of  the  empirical 
formula  given  in  the  original  text,  the  following  exact  formula 
is  obtained: 

ft'»s.nr  =  *V  ,.  +  {?<*»,  o-fr'.u  »)  «A-4I 

A  further  refinement  of  the  present  analysis  is  the  discovery 
of  an  exact  integral  for  the  cavitation  boundary,  without 
invoking  the  “mean  transport  approximation”  1 1 1  yielding 


fi,-  <L(’V  dr,  (V?» 


where  the  subscript  **B"  refers  to  the  cavitation  boundary 
either  in  the  break-up  or  in  the  Fill-back  region.  For  the  break 
up  boundary,  the  right-hand  side  of  equation  (A-5)  reduces  to 
Q,  fQh  in  agreement  with  equation  ( 17)  by  the  limiting  process 
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Vn  —  n/>-  For  the  fill-back  boundary,  the  subscript  “fl”  is 
replaced  by  *T*  for  (f,  H,  if >.  ( Hh ,  if* I  refer  to  the  break-up 
boundary  at  Equation  <A:5)  thus  becomes  a  tran¬ 
scendental  relationship  between  f,  and  and  can  be  solved 
only  numerical  y.  The  actual  discrepancy  between  equation 
(21),  which  is  based  on  the  “mean  transport  approximation,” 
and  equation  (A-5)  is  quite  minute  in  all  physically  significant 
aspects. 

It  is  understandable  that  the  analytical  content  of  the 
present  work  should  be  regarded  as  being  too  complex  for 
practical  use.  However,  it  is  because  of  the  fully  analytical 
treatment  that  very  complete  numerical  results  can  be 
compiled.  In  view  of  their  universal  applicability,  tabulat.ons 
for  the  rull  range  of  U ,  *  fare  furnished  with  this  closure.  It  is 
hoped  that  the  availat  ty  of  the  numerical  tables  will  en¬ 
courage  judicial  use  of  these  results  in  pertinent  engineering 
activities.  The  headings  used  in  these  tables  are  defined  as 
follows: 
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The  interesting  experimental  results  reported  by  Dr.  Etsion 
w  ill  no  doubt  enhance  our  knowledge  of  the  cavitation  process 
in  fluid  film  bearings.  This  author  looks  forward  to  the  full 
publication  of  his  findings.  In  the  context  of  the  present 
discussion,  this  author  would  suggest  that  hydrodynamic 
action  in  the  gaseous  phase,  along  with  inertia  effects  and 
mass  transfer  related  to  gas  phase  solubility  may  be  all 
associated  with  the  interesting  observations  reported  by  Dr 
Etsion.  It  should  be  clear  that  the  details  of  the  cavitation 
process  must  depend  on  other  thermodynamic  and  fluid 
dynamic  parameters  which  are  overlooked  in  conventicanl 
lubrication  studies.  For  this  reason,  in  terms  of  details  of  the 
cavitation  process*  universality  among  journal  hearings, 
seals,  and  dampers  will  cease  to  exist. 
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Dynamic  Analysis  of  Rupture  in 
Thin  Fluid  Films.  I— A  Noninertial 
Theory' 

A  noninertial  theoretical  model  for  the  dynamics  of  film  rupture  has  been  for 
undated.  Under  the  transient  condition,  movement  of  the  rupture  boundary  h 
governed  hv  the  condition  of  flow  continuity  between  the  film  flux  and  adhered  film 
transport  in  the  cavitation  domain.  The  traditional  Swift -Srieber  condition  for  film 
breakup  is  shown  to  be  valid  upon  reaching  steady-state.  Generalization  is  extended 
to  allow  consideration  of  two  sliding  surfaces  and  the  pure  squeeze-film.  The 
possibility  of  subcavitation  film  pressure  is  shown  to  result  in  dry  regions  in  the 
cavitation  domain. 


1  Introduction 

In  the  lubricant  film  of  a  self-acting  bearing,  surface  sliding 
in  the  presence  of  spatial  variation  of  the  bearing  gap  con¬ 
stitutes  the  hydrodynamic  wedge  action,  which  is  responsible 
for  the  load  carrying  pressure  in  the  bearing  film.  If  surface 
sliding  is  in  ihc  direction  ot  gap  divergence,  (hen  the  pressure 
in  the  lubricant  film  tends  to  decrease.  With  the  proper  gap 
distribution,  the  reduced  pressure  in  the  diverging  gap  can 
combine  with  the  elevated  pressure  in  the  convergent  gap  to 
form  the  overall  load  capacity  of  the  bearing  as  illustrated  by 
the  classical  Sommerfcld  solution  of  the  infinitely  long 
journal  bearing  |l|.  Although  a  reduced  film  pressure  in  a 
bearing  may  not  be  undesirable,  it  hardly  ever  happens  in 
reality.  It  is  usually  found  that  the  pressure  in  a  conventional 
lubricant  film  cannot  fall  substantially  below  the  atmospheric 
level.  Presumably,  lubricant  used  in  real  environments  always 
contains  some  dissolved  ambient  gas,  volatile  constituents  of 
moderately  high  vapor  pressure,  and  sometimes  even  en¬ 
trained  gas  phase  in  suspension.  Thus,  as  surface  sliding 
enters  into  a  divergent  region,  the  void  fraction  in  the  bearing 
gap  will  rapidly  increase  through  a  combination  of  physical 
and  chcm  cal  processes.  As  a  result,  this  region  of  high  void 
fraction  will  prevent  further  reduction  of  the  film  pressure 
according  to  the  hydrodynamic  lubrication  theory,  instead, 
expansibility  and  a  relatively  low  viscosity  of  the  gas  phase 
allow  a  zone  of  nearly  uniform  pressure  to  be  formed.  Among 
workers  in  the  field  of  lubrication  research,  this  process  is 
interchangeably  known  as  cavitation  or  film  rupture. 

A  totally  rigorous  treatment  of  the  film  rupture  process 
should  give  some  attention  to  the  thermodynamic  and 
chemical  behaviors  of  the  liquid  lubricant.  The  present  work 
is  a  less  ambitious  undertaking  by  assuming  a  purely 
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mechanistic  viewpoint.  Allowing  sufficient  time  lapse,  phase 
equilibrium  may  be  considered  to  have  been  established.  The 
pressure  level  in  the  cavitated  lubricant  film  should  ir  some 
way  be  related  to  the  physical  chemistry  of  the  lubricant,  and 
may  be  regarded  as  a  constant  parameter  in  the  context  of  a 
mechanistic  study.  The  detail  flow  structure  of  the  ruptured 
film,  viewed  in  the  scale  of  the  film  thickness,  is  quite  in 
tricate.  Liquid-solid-gas  interfacial  effects,  free  surface 
curvature,  and  fluid  inertia  are  simultaneously  interacting 
with  the  viscous  stress  which  is  dominating  in  the  remote 
regions  (from  the  rupture  boundary).  Meniscus  feature** 
and/or  striation  fingers  are  merely  particular  features  of  the 
flow  field  here.  The  present  work  shuns  involvement  in  thb 
fine-scale  point  of  view  and  will  deal  mainly  with  gros^ 
features  on  either  side  of  the  rupture  boundary.  Surface 
tension,  for  a  film  thickness  of  the  submicron  level,  can  also 
manifest  itself  in  a  subcavitation  pressure  in  the  immediate 
upstream  of  the  rupture  boundary.  In  the  gross  scafe,  one  can 
accommodate  this  phenomenon  by  assigning  a  pressure  jump 
at  the  rupture  boundary.  For  common  situations  associated 
with  bearings  and  dampers,  the  amount  of  this  pressure  jump 
is  quite  small  when  compared  with  the  film  pressure 
associated  with  hydrodynamic  wedge  and/or  squeeze  effects, 
and  thus  can  justifiably  be  neglected.  Pressure  level  alone, 
however,  does  not  uniquely  define  the  theoretical  problem  of 
a  ruptured  lubricant  film;  thus,  the  zero  pressure  gradient 
condition  (2,3)  has  been  regarded  as  the  rigorous  rupture 
condition.  Use  of  the  latter  condition  to  analyze  a  one 
dimensional  bearing  film  (infinitely  tong)  is  a .  common 
practice.  Two  dimensional  problems  are  less  frequently 
treated  with  the  same  degree  of  care,  due  to  the  extreme 
complexity  in  computational  undertaking  required.  The  work 
of  Jakobsson  and  Floberg  (4)  is  an  exemplifying  exception;  it 
brought  out  the  need  to  satisfy  the  continuity  condition 
throughout  the  cavitation  domain.  In  lieu  of  satisfying  the 
rigorous  rupture  condition,  it  is  a  common  approximation  to 
treat  the  film  rupture  problem  with  the  Gtimbcl  hypotheses 
|5|;  such  that  ihc  film  pressure  is  first  solved  with  no  regard 
,  for  film  rupture,  then  all  subambient  pressure  values  are 
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replaced  by  the  ambient  pressure.  Such  an  approach  is  par¬ 
ticularly  popular  in  connection  with  the  short  bearing  ap¬ 
proximation  (6)  under  the  dynamic  condition.  More  recently, 
♦he  Swift-Stieber  condition  of  film  rupture  has  been  adapted 
to  study  short  bearing  configurations  [7-IOJ;  in  comparison 
with  the  Gumbel  approximation,  significantly  different 
results  were  obtained  (9,101. 

While  there  is  now  a  thorough  understanding  of  the  issues 
associated  with  a  steady -state  film  rupture  process,  relatively 
little  attention  has  been  given  to  the  time-dependem  coun¬ 
terpart.  E  rod  |11)  proposed  a  computational  procedure  for 
the  time-dependent  cavitation  according  to  a  derivation  due  to 
Olsson  [12].  Fill  fraction  of  the  gap  was  identified  as  the  state 
variable  in  a  ruptured  film,  and  any  movement  of  the 
boundary  between  the  full-film  and  the  ruptured  film  was 
treated  as  a  time  rate  of  change  of  the  fill  fraction.  Choosing 
to  rely  on  a  single  algorithm  for  both  the  full  film  and  the 
cavitation  domains,  the  concept  of  fill  fraction  was  retained 
for  the  full  film  via  the  use  of  fluid  compressibility  so  that 
pressure  elevation  above  the  cavitation  pressure  was 
associated  with  a  fill  fraction  in  excess  of  unity.  In  the  present 
work,  the  problem  of  a  lime  dependent  film  rupture  problem 
is  examined  analytically  in  a  manner  consistent  with  the 
established  lubrication  theory  of  the  full  film. 

Primary  emphasis  of  the  present  work  will  deal  with  the 
self-acting  bearing,  in  which  a  closed  cavitation  domain  limits 
the  lower  bound  of  the  film  pressure.  Wedge  action  and 
squeeze  film  together  determine  the  film  pressure  according  to 
the  lubrication  theory.  Tne  sliding  motion,  traversing  between 
rupture  boundaries,  dominates  the  flow  distribution  within 
the  cavitation  domain. 

Film  rupture  is  known  to  exist  in  squeeze  bearings  with  or 
without  simultaneous  sliding.  In  particular,  a  complete 
mathematical  equivalence  exists  between  a  squeeze-film 
damper  in  a  circular  whirl  condition  and  a  statically  loaded 
journal  bearing,  including  the  global  structure  of  cavitation. 
The  absence  of  a  sliding  motion,  however,  brings  about  a 
different  situation  in  the  flow  balance  at  the  rupture  boun¬ 
dary.  Therefore,  film  rupture  dynamics  in  the  squeeze  film 
will  be  given  separate  attention. 


Since  film  rupture  is  generally  regarded  as  a  mechanism  to 
limit  an  excessive  reduction  of  film  pressure,  the  cavitation 
pressure  is  usually  assumed  to  be  lower  than  the  film  pressure 
in  the  interior.  However,  there  is  no  fundamental  argument  to 
rule  out  a  subeaviration  film  pressure  in  the  absolute  sense, 
particularly -under  a  dynamic  environment.  Therefore*  in  the 
present  study,  the  consequence  of  a  subcavitatior  film 
pressure  will  be  examined. 

As  the  first  of  several  articles,  the  author  intends  to  ex 
pound  problems  related  to  the  dynamics  of  ruptured  fluid 
films;  he  chooses  not  to  invoke  computed  example^.  It  i*. 
hoped  that  broadly  applicable  conclusions  can  be  firrnK 
established  in  the  qualitative  sense.  Subsequent  efforts  can 
then  be  properly  channeled  to  seek  specific  results  of 
significance. 

Five  major  assumptions  will  be  observed  in  subsequen 
derivations.  They  arc  briefly  stated  as  follows: 

•  Standard  approximations  of  the  lubrication  theory  arc 
retained  for  the  unruptured  domain. 

•  A  constant  and  uniform  pressure  exists  in  the  ruptured 
domain. 

•  A  homogenous  flow  structure  is  assumed  in  the  ruptured 
domain. 

•  An  abrupt  change  in  the  flow  structure  occurs  at  the 
rupture  boundary.  Surface  tension  is  neglected,  and  pressure 
continuity  is  assumed. 

•  Inertial  effects  are  neglected  in  the  ruptured  domain. 

2  The  Self-Acting  Bearing 

The  self-acting  bearing  of  the  conventional  type  is  con¬ 
sidered  here.  The  adjective,  “conventional/*  is  emphasized  to 
exclude  the  absence  of  a  surface  velocity  and  such  special 
bearings  which  have  nonvanishing  sliding  velocities  at  both 
surfaces. 

With  only  one  sliding  surface,  the  flow  field  in  iht 
cavitation  domain  will  be  modeled  as  a  layer  of  fluid  which  is 
adhered  to  the  sliding  surface.  The  adhered  fluid  layer  will  be 
called  the  “adhered  film/*  No  interaction  between  tin- 
adhered  film  and  the  stationary  surface  is  allowed  tc  cause 


Nomenclature 


B  ~  a  functional  form  of  the 
external  boundary  condition 
C  =  nominal  gap 
e  ~  journal  bearing  eccentricity 
h  =  gap  or  film  thickness 
ha  =*  thickness  of  adhered  film 
hf  =  film  thickness  at  a  rupture 
point 

nf  =  unit  outward  normal  of  the 
rupture  boundary 
p  =  film  pressure 
pc  =  pressure  in  the  cavitation 
domain 

p,  ~  lubricant  supply  pressure 
Q  «  film  flux  . 

rB  =  vector  coordinate  of  the 
external  boundary 
rr  «  vector  coordinate  of  the 
rupture  boundary 
R  »  journal  radius 
f  =  lime 

th  a  birth  time  of  adhered  film 


tf  »  final  time  of  adhered  film 
tQ  ~  characteristic  time 
U,  =  sweeping  speed  of  the  rupture 
boundary 

Ur  =  sweeping  velocity  of  a 
rupture  point 

V  -  sliding  speed 

V  *  sliding  velocity 

x  =  circumferential  coordinate 
on  a  journal  bearing  surface 
x0  =  characteristic  circumferential 
coordinate 

.v*  =  circumferential  coordinate  of 
a  rupture  point;  function  of 
(z,t) 

z  -  axial  coordinate  of  a  journal 
#  bearing 

z,  =  axial  coordinate  of  a  rupture 
point;  function  of  (x,  /) 

Symbols 

0  «  inclination  of  the  rupture 
boundary  to  the  surface 
velocity 


n  =  fluid  viscosity 
a  =  striation  fraction 
u  =  rotation  speed 

Subscripts 

a  =  adhered  film 
b  =  breakup  (of  full-film) 
c  =  cavitation 
/  =  fill-back  (of  full-film) 

/  =  rupture  initiation 
0  =  reference  point  of  the 
transport  characteristics 
p  =  related  to  the  pressure 
(Poiseuille)  effect 
r  *  rupture  boundary 
R  =  film  reestablishment 

V  =  related  to  the  sliding  effect 
W  =  widest  point  of  cavitation 

V  *  two-dimensional  gradient 

operator  on  the  bearing 
surface 

(  *  )  =  variable  (  )  view  in  a  frame 
of  reference  fixed  to  sut  face  1 
1.2  =  surfaces  I  and  2 


•When  Mit-uript  ,V*  i\  replaced  by  or  a  breakup  point  or  a  fill- 
back  point  is  ndicated. 
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Fig.  1  Gaomatry  of  an  ancloaad  cavitation  domain 

pressure  variation  in  the  cavitation  domain.  As  a  con¬ 
sequence,  upon  neglecting  any  inertia  effects,  the  adhered 
film  can  he  regarded  to  he  fully  quiescent  relative  to  the 
sliding  suracc. 

This  model  can  be  generalized  to  consider  the  bearing  with 
two  sliding  surfaces  provided 

•  a  franc  of  reference  <an  be  used  so  that  one  of  the  two 
surfaces  is  not  sliding  with  respect  to  the  moving  frame  of 
reference, 

•  the  second  surface  has  a  nonvanishing  sliding  velocity 
relative  to  the  moving  frame  of  reference,  and 

•  the  sum  of  the  sliding  velocities  does  not  vanish. 

Rules  to  be  followed  in  order  to  achieve  this  generali/at  ton 
will  be  disci  ssed  at  a  later  time.  The  immediate  derivations 
will  apply  only  to  the  case  of  a  conventional  self-acting 
bearing  witli  a  single  sliding  surface.  The  presence  of  (time 
dependent!  squeeze  effect,  in  addition  to  the  wedge  effect  of 
the  sliding  surface,  however,  may  be  considered  with  the 
usual  inclusion  of  the  squeeze  term  in  the  treatment  of  the 
film  pressure  and  by  usii  g  the  instantaneous  gap  for  the  fluid 
film  at  the  rupture  boundary. 

(A)  Geometrical  Description  of  a  Closed  Cavitation 

Domain.  If  the  cavitation  domain  should  be  enclosed  by  a 
full-film  on  alt  sides,  as  in  the  case  of  a  journal  bearing  with 
both  ends  flooded  or  supplied  with  fluid,  the  rupture 
boundary  can  be  divided  into  upper  and  lower  branches  for 
the  corner iencc  of  analvtica!  treatment.  In  the  case  of  a  self¬ 
acting  bearing,  each  branch,  in  turn,  can  be  divided  into  a 
breakup  portion  and  a  fill-back  portion.  The  four  parts  of  the 
rupture  boundary  can  be  identified  with  the  inclination  of  its 
outward  normal  (pointing  into  the  cavitation  domain)  relative 
to  the  surtax c  velocity.  According  to  the  present  convention, 
the  surface  velocity  being  aimed  along  the  positive  v-axis,  the 
inclination  of  the  outward  normal  has  the  following  indicated 
property  on  each  of  the  four  parts  of  the  rupture  boundary: 

Breakup  Till-back 

Upper  branch  -acute  obtuse 

Lower  branch  +  acute  +  obtuse 


As  illustrated  in  Fig.  I ,  the  upper  and  lower  branches  of  the 
rupture  boundary  part  from  each  other  at  the  rupture 
initiation  point  *7  "  and  continue  on  their  respective  breakup 
portions  moving  laterally  away  from  each  other  Each 
breakup  boundary  becomes  more  and  more  aligned  with  the 
surface  velocity  vector  until  it  reaches  point  •‘MV*  where  the 
boundary  grazes  the  surface  velocity  vector.  Beyond  point 
"W"  the  boundary  turns  back  in  becoming  a  fill-back  line. 
The  fill -back  lines  eventually  meet  again  at  the  film 
reestablishment  point  "W.” 


While  a  rectangular  pi.  orm  is  illustrated  in  Fig.  I.  these 
features  of  the  cavitation  domain  and  its  bound: v  are 
generally  valid.  For  example,  the  bearing  plan  form  nia>  be 
annular,  such  that  the  s  ir  face  velocity  is  a  curved  field. 

(B)  Fluid"  Transport  in  the  Cavitation  Domain,  in  a 

cavitatcd  or  ruptured  domain  of  a  self-acting  hearing,  the 
concept  of  an  adhered  film  may  be  used  to  describe  th.‘  fluid 
transport  process.  Explicitly,  the  adhered  film  thickness  ma> 
be  designated  h,\  which  can  also  be  used  for  a  striated  film 
and  be  related  to  the  vl  nation  fraction  u  according  to 

ht.  =  oh  i\) 

Thus,  subsequent  discussions  regarding  fluid  transport  are 
applicable  to  both  an  adhered  film  in  reality  or  the  striated 
film  by  analogy.’ 

One  may  reasonably  postulate  that  the  gaseous  const  iuicnp 
within  the  void  regions  do  not  influence  the  motion  of  the 
adhered  film.  Accordingly,  the  fluid  transport  law  for  the 
adhered  film  can  be  stated  as 

This  is  a  first  order  partial  differential  equation  of  thi 
hyperbolic  type,  which  possesses  one  real  characteristic 
equation  [13).  For  the  particular  situation  on  hand,  us  ng  the 
uniform  field  V  of  the  cylindrical  journal  bearing  fot 
illustration,  the  characteristic  equation  is 

v -,v0  =  T(/w„)  (3 

The  axial  coordinate  does  not  appear  explicitly  in  *hi 
characteristic  equation  because  equation  (2)  makes  nc 
reference  to  the  direction  transverse  to  the  surface  vt  loots 
For  the  same  reason,  the  characteristic  equation  is  nicauingfu 
only  when  the  axial  coordinate  is  fixed.  may  be  zero,  ihei 
a„  refers  to  an  element  of  the  adhered  film  in  its  initial  field  h 
(.v„,  C,  0).  Or.  for  t„  -  th> 0.  v„  is  to  be  identified  with  the 
breakup  point  xh{z*  M  on  the  rupture  botndan .  The 
characteristic  equation  is  the  trajectory  of  an  element  of  :hc 
adhered  film  in  the  \  -(  plane  for  a  fixed  axial  location.  Thi 
trajectory  terminates  as  the  element  of  the  adhered  filir  meet 
a  fill-back  point  on  the  rupture  boundary  .v,  (c.  M. 

(C)  Sweeping  Motion  of  the  Rupture  Bowtdan .  The 

rupture  boundary.  whi;h  separates  the  full-film  and  the 
cavitation  domains,  may  not  be  stationary.  Its  mo  ion  i 
dictated  by  the  balance  of  flux  between  the  full-film  and  the 
adhered  film.  By  definition,  the  sweeping  velocity  l\  is  aimcc 
along  the  outward  normal  n,.  As  illustrated  in  Fig.  2,  which 
depicts  the  displacement  by  the  sweeping  motion  for  an  in 
finitesimal  duration  A/,  continuity  consideration  yields 

*/,<*, )<▼#».+  ~  V -*„%■).«,  (4) 

The  first  term  on  the  right-hand  side  is  the  Poiscuillc  flow 
Because  the  rupture  boundary  is  an  isobar,  it  is  always  per 
pcndicular  to  the  Poiscuillc  flow.  If  the  film  pressure  is  highei 
than  the  cavitation  pressure,  then  the  Poiseuillc  flow  tends  to 
cause  the  rupture  boundary  to  contract;  the  reverse  is  true  it 
the  film  pressure  is  lower.  The  second  and  third  terms  are. 
respectively,  the  Couettc  flow  and  the  transport  of  the 
adhered  film.  Being  directionally  controlled  by  the  surface 
velocity,  they  contribute  to  the  motion  of  the  r  ipture 
boundary  through  the  normal  component  (to  ihc  r  jpttirt 
boundary)  of  the  surface  velocity.  Depending  on  the  balance 


1  A  fundamental  tliainction  bciwcen  an  adhere  J  film  and  a  Mr  au\l  f  ton  w  ih 
ability  of  the  latter  to  sustain  a  shear  stress.  This  distinction  is  important  onl\ 
when  attempts  to  calculate  friciion  arc  made,  but  is  irrelevant  u»  <Mti.it ten,  o* 
prcsviirc.  How,  and  the  ncotncirv  of  the  rnpfure  botmdar> . 
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Fig.  2  Di*pl*cem*nt  ol  a  twtaping  rupture  point 


between  the  Couettc  flow  and  the  adhered  film  transport,  as 
well  as  the  inclination  of  the  surface  velocity  relative  to  the 
rupture  boundary,  the  net  effect  of  these  two  terms  may  also 
cause  the  rupture  boundary  to  contract  or  to  enlarge. 

For  a  cy  indrica)  journal  bearing,  the  inclination  of  the 
rupture  boundary  is  0  =  t/2  -  tan  1  (Ax ,/Az)  and  that  of  the 
outward  normal  is  (3  *  t/2,  the  alternative  signs  of  the  latter 
angle  apply,  respectively,  to  the  upper  and  lower  branches  of 
the  rupture  boundary.  Hence  equation  (4)  becomes 

*(-(£)©.  <» 
The  foregoing  equatior  encompasses  the  special  case  of  a 
stationary  tupture  boundary,  for  which  the  breakup  point 
requires 


(S).-(S).-  — 

while  the  fill-back  point  requires 


Equation  (6b)  may  be  recognized  as  the  Swift-Stiebcr 
condition  |2,3J.  The  steady-state  problem  for  a  journal 
bearing  of  finite  length  was  solved  by  Jakobsson  and  Floberg 
[4).  The  short -bearing  version  of  the  steady-state  problem  was 
treated  by  Wakurai,  et  al.  (7),  with  respect  to  the  flow 
requirement  of  a  short  journal  bearing,  by  Findlay  (8)  con¬ 
cerning  the  opera  >n  of  a  wavy  face  seal,  and  by  Pan  (9,  10] 
regarding  the  load  characteristics  of  a  short -bearing  or  a 
squeeze-film  damper. 

For  a  time -dependent  problem,  an  important  issue  is  the 
appropriate  relationship  between  A*,  the  film  thickness  of  the 
full-film  approaching  the  breakup  boundary,  and 
the  thickness  of  the  adhered  film  immediately  downstream. 
This  question  must  be  resolved  in  the  context  of  providing  a 
consistent  link  between  the  full-film  and  the  adhered  film 
across  the  breakup  boundary.  Since  the  adhered  film  in¬ 
troduces  the  point  of  view  of  a  characteristic  line  at  a  fixed  z, 
equation  (3\  t  is  useful  to  rewrite  equation  (5)  into 


■((£)©,  (5) 


The  subscript  “r”  indicates  that  equation  (8)  is  in  a  form 
which  is  applicable  to  cither  a  breakup  point  or  a  till  back 
point.  However,  for  the  breakup  point,  a  relationship  between 
hh  and  h(l  (zM  th)  must  somehow  be  established;  while  for  the 
fill-back  point,  ht,  (.v,.  /,)  is  lo  be  consistent  with  the  traits 
port  law  for  the  adhered  film. 

Actual  implementation  of  equation  (8)  would  become 
difficult  in  the  neighborhood  of  point  W\  where  (A.\r/Ac) 
becomes  unbounded.  This  difficulty  can  be  avoided  by 
considering  the  rupture  boundary  as  z,  (x,  t)  in  this  vicinity. 
Then,  instead  of  equation  (8),  one  can  w  rite 

dzr 

=  T  Ur  seed 
dt 

♦[(£)(£). -"(1 -*>](*■)>  «* 

(D)  World  Diagrams  of  the  Adhered  Him.  The  transport 
process  of  the  adhered  film,  as  described  by  the  characteristic 
equation,  can  be  effectively  presented  in  terms  of  the  world 
diagram  which  traces  the  space- time  relationships  of  pertinent 
phenomena.  As  an  illustrative  example,  a  model  problem  will 
be  described  to  depict  the  evolution  of  the  cavitation  domain 
from  that  of  a  “half  Sommerfeld  film’*  to  the  steady  state 
solution  which  rigorously  satisfy  the  Swift-Stieber  condition 
and  the  adhered  film  transport  jaw.  For  this  purpose,  the 
bearing  gap  h  (x,  y)  is  assumed  to  be  time-independent.  The 
initial  state  of  the  full-film  domain  is  that  portion  of  the  “full 
Sommerfeld”  solution  where  p  (x,  v)  2  p ,  being  the 

assumed  cavitation  pressure.  Within  the  rupture  boundary  x, 
(z,  t  =  0)  which  encloses  the  subcavitation  domain  of  the  lull 
Sommerfeld  solution,  a  suitable  initial  field  of  the  adhered 
film  is  presumably  known,  i.e.,  ha  (x,  z,  /  *  0)  has  been 
specified.  For  a  fixed  value  of  z,  one  is  thus  given  the  com¬ 
plete  distribution  of  ha  (/  *0)  over  the  span  xh  (z,  f  =  0)  <  x 
s  Xj  (z,  (  -  0).  A  series  of  characteristic  lines  can  be  drawn 
in  th  ex  -  f  plane  as  described  by  the  characteristic  equation 

x(/)  -x(0)=  Vt  (10) 

Along  each  characteristic  line,  ha  retains  its  value  us 
specified  in  (x,  z,  t  =  0).  The  world  diagram  is  bounded  on 
the  left  by  the  history  of  the  breakup  point  x„  (z,  f)  and  on  the 
right  by  the  history  of  the  fill-back  point  xf  (z,  f).  xh  and  xf 
are  obtained  by  (he  time-domain  integrations  of  (dxh/di)  and 
(dXf/dt)%  respectively. 

The  left  boundary  of  the  world  diagram,  jc„,  immediately 
brings  up  the  question  regarding  the  proper  value  of  hQ  (x*.  z, 
th),  in  order  to  initiate  the  analysis  of  the  transport  process. 
Previously,  it  was  shown  that  a  stationary  breakup  boundary 
requires  that  both  equations  (6a)  and  (6b)  be  satisfied. 
Violation  of  equation  (6b)  causes  the  breakup  boundary  to 
move.  However,  so  long  as  the  breakup  boundary  lags  behind 
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i he  sliding  surface  (subcritical  feeding),  equation  (6/?)  can  be 
let  allied;  l  bus, 


h „  (xh,z*ih) 


for  subcritical  feeding 


(lltf,  b) 


ill  V  6^x  /  L  V  dz  '  h  '  &Z  '  V  Av  '  **  J 

This  sci  of  results  cannot  be  universally  valid,  because,  in 
the  vicinity  of  point  W  where  dxh/dz  —  oo#  equation  (l  ib) 
wouid  become  unbounded.  Before  the  latter  condition 
prevails,  \xt,/df  would  exceed  V  so  that  whatever  may  be  h „ 
t  v...  c,  r„),  the  motion  of  the  breakup  point  would  sweep 
ahead  the  adhered  film  such  that  ha  would  have  to  be 
associated  with  the  characteristic  field  of  prior  time,  i.e.,  the 
adhered  film  at  the  point  immediately  downstream  of  the 
breakup  point  has  existed  prior  to  the  approach  of  the  rupture 
boundary.  Thus  the  last  term  on  the  right-hand  side  of 
equation  (8),  which  was  omitted  on  account  of  equation  (I  I  u), 
would  have  to  be  retained.  Thus,  the  subcritical  feeding 
condition  must  obey  the  constraint 


(£)[(£). (£MS).H  «» 

A  super  critical  feeding  condition  is  indicated  if 


OlOAtnt)  I-— •>“ 


(1J> 


Substituting  into  equation  (8)  with  hd  identified  with  the 
prior  time  characteristics,  one  obtains 


dV/.  _  y 

<U  2  ih„-hu) 


(14) 


whici  assures  that  the  breakup  point  would  indeed  sweep 
ahead  the  adhered  film  regardless  of  the  specific  value  of  hu. 
Since  the  condition  of  supercritical  feeding  is  most  likely 
induced  by  a  large  magnitude  of  (d.vr/dc),  sweeping 
moti  m  should  be  calculated  in  terms  of  (dzr/9t)  according  to 
cqua  ion  i9).  However,  the  test  for  the  critical  feeding  con- 
dmoi  should  still  be  performed  in  terms  of  equation  (13). 
Subsequently, 


a;. 

it 


-  (t->[(s>. -£).(s)l. 


for  subcritical  feeding 


— *■>  ■(-(£)©. 

•[(«)(«). -Kt-OKt)) 


for  supercritical  feeding 

(1 5fl,b) 


The  right  boundary  of  the  world  diagram  has  a  motion 
described  by  equation  (8)  if  one  merely  replaces  the  subscript  r 
by  /.  It  may  be  regarded  to  be  the  sum  of  a  pressure  com¬ 
ponent  idxf/dt)r  and  a  sliding  component  [9xf/Bt) , : 

«*V,  (9x,\  (*X'\ 

it i  "  V  Hi  K  V  9,  /, 


(0)  Super -Critic ot  FM#ng 


! 


(Cl  For  motion  of 

StMtfy- Stote  Watt* 


Fig.  3  World  diagram*  of  r upturn  transition 


/<*'•,  \  _ _ v _ ((d/,\  (3h )  - \ 

V  ,V K  OS-/,  '  J:' 

(Av')  =  I-.--*'-  )  (Ito.ft.c) 

V  fit  '«  2 {hf-hj  V  2 ihf-hj' 

If  ihc  film  pressure  is  assumed  to  be  above/?  at  all  times, 
[hen  (dx,/dt)r  is  always  negative.  Now  that  {dxf/dt)t  is 
always  less  than  V\  therefore  (dx,/i )/)  is  always  i Igebraically 
less  titan  V  so  that  one  can  be  certain  that  the  fill -tack  point  is 
always  overtaken  by  the  adhered  film.  Again,  where  the 
magnitude  of  \9xffbz)  is  large,  it  is  better  to  treat  he  fill-back 
boundary  as  z,  (x,  r)  in  numerical  implementation. 


iV 


(17) 


Representative  world  diagrams  showing  the  transition  of 
the  cavitation  domain  in  a  journal  bearing  film  from  the  half 
Sommcrfcld  state  to  the  steady-state  Swift-Stieber  state  are 
given  in  Fig.  3.  Each  world  diagram  describes  the  history  of 
the  adhered  film  at  a  particular  axial  location.  At  the  bottom 
of  Fig.  3.  solid  and  dashed  outlines,  tespectively,  depict  the 
initial  and  final  rupture  boundaries.  Three  horizortal  lines  are 
drawn  through  both  boundaries.  zA  is  near  the  central  region 
of  the  cavitation  domain  In  this  case,  the  breakup  point  feeds 
at  a  subcritical  condition.  The  fill-back  point  moves  under  the 
influence  of  the  initial  field  ha  (* ,e4,  0)  until  it  arrives  at 
point  ,  where  the  earliest  characteristic  line  Irom  xh  in¬ 
tersects  xf .  Both  Jr*  and  xf  have  xertical  asymptotes  which 
mark  the  steady-state  boundary  points.  zB  is  dose  to  the 
initial  zero  slope  point.  The  breakup  point  moves  ahead  of  its 
own  characteristic  lines.  The  fill -back  point,  xvith  bxf/bt 
dominated  by  its  pressure  component*  moves  backward  and 
meets  the  breakup  point  at  *•£.”  The  world  diagi  a  n  ends  here 
because  the  rupture  boundary  has  swept  below  c, .  zr  traces 
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toim  W  of  the  final  rupture  boundary.  The  breakup  point 
now*  forward  rapidly  at  first  *  albeit  at  a  suberitical  feeding 
ondition.  The  fill-back  point  moves  back  quite  rapidly  in  the 
uitial  period,  dominated  bv  Ihe  pressure  component  under 
he  influence  of  the  initial  field  h „  (a,  z{  ,  0)  until  point  *'S,  ," 
'.here  a,  meets  the  earliest  characteristic  line  from  xh.  In  the 
.tic  period,  as  both  x„  and  a*,  approach  the  final  location  of 
Hfim  M ,  the  magnitude  of  the  normal  pressure  gradient  in  the 
luid  film  at  the  rupture  boundary  reduces  to  a  negligible 
tvel:  they  both  turn  toward  the  vertical  direction,  reflecting 
lowing  down  of  the  sweeping  motion  and  eventually  they 
.oproach  the  same  asymptote  and  coalesce  into  the  final 
ocat  ion  of  point  W. 


•  Calculate  the  sweeping  rates  of  the  rupture  boundary, 
using  equations  (8)  and  (9)  as  may  be  appropriate.  I  nitially,  a 
starling  di  cribution  of  the  adhered  film  thickness,  ha  (r, 
f  =  0)  in  the  cavitation  domain  is  to  be  furnished,  e.g., 
equation  (18). 

•  Step  in  time  to  update  the  geometry  of  the  rupture 
boundary.  Subsequently,  the  full  cycle  of  computations  can 
be  repealed  as  many  limes  as  necessary.  To  enhance  the 
ovetall  accuracy  of  the  time-stepping  operation,  the  Runge- 
Kutia  method  and/or  the  predictor  corrector  method  1 14] 
may  be  adopted.  The  size  of  time  step  should  reflect  the 
resolution  required  to  depict  the  time-dependence  of  o. 


(K)  The  Initial  Kidd.  The  true  initial  period  of  void 
duration  and  growth  must  necessarily  involve  mass  transfer 
’Ctween  liquid  and  gas  phases  so  that  thermodynamic  and 
hemicnl  parameters  may  not  be  negligible.  The  purely 
nechamsfic  model,  therefore,  may  not  be  applicable  until  a 
.tier  time,  when  the  status  quo  of  a  relatively  extensive  void 
taction  prevails  in  the  cavitation  domain.  The  initial  field  of 
ic  thickness  distribution  of  the  adhered  film,  hu  ( v,  c.  0), 
uust  be  reconcilable  with  the  initial  film  pressure.  If  one 
u’cepts  the  constraint  that  the  film  pressure  cannot  fall  below 
ie  cu*  it  at  ion  level,  then  it  is  necessary  to  exclude  the 
ossibi.ity  of  a  fully  filled  gap  as  the  initial  condition,  which 
eqjires  h,  (a\  c.  0)  -  h  ( v,  c,  f  —  0).  For  this  reason,  the 
pseudo  steady -state  condition  to  be  defined  as  follows  is 
.  ropos.'d  to  be  the  initial  state  of  the  adhered  film.  Given  the 
description  of  the  rupture  boundary  at  the  initial  instant.  x, 

( z ,  0),  he  pseudo  steady-state  thickness  of  the  adhered  film  is 

/i.  (v,. <v<\,.:,0)  -  (18) 

It  should  be  recognized,  however,  the  constraint  that  the 
Mini  picssuie  cannot  be  subambient  may  not  be  realistic. 

I  here  na\  be  circumstances  tor  which  such  an  assumption 
should  be  relaxed.  This  issue  will  he  further  discussed  at  a 
'  .iter  section. 

(F )  Interfacing  with  Film  Pressure  C  alculation.  Analysis  of 
he  adhered -film  transport  problem  cannot  be  performed 
without  dealing  with  the  surrounding  fluid  film  pressure  at  the 
•aroc  i*me.  In  fact,  the  starting  point  is  most  likely  a  half 
''OinnuTicld  type  pressure  field  which  is  the  solution  of  the 
.oinentional  Reynolds'  equation.  Steps  of  computation 
required  to  perform  the  overall  analysis  are  outlined  in  the 
following. 

•  Solve  the  appropriate  form  of  Reynolds’  equation,  e.g., 

V*l-(|2„)V/,+  2V  *)“°  <l9) 

Requisite  boundary  conditions  include  both  external  ones, 
e.g,.  tfosc  associated  with  hearing  ends  and  feed  slots,  which 
car  he  expressed  in  the  form 

«|r„;  />./?,  |=0  (20) 

rH  is  the  vector  coordinate  of  the  external  boundary  and  />,  is 
(he  applicable  supply  pressure,  and  also  the  pressure  of  the 
tupture  boundary 

/Mr,,  /)=/>,  (21) 

r,  is  the  vector  coordinate  of  the  rupture  boundary.  An  initial 
estimate  of  ihe  rupture  boundary  may  be  determined  from  the 
half  Sommerfeld  solution.  Pressure  gradient  at  the  rupture 
boundary,  Vp  (r,,  /),  is  to  be  established. 

•  Subdivide  the  rupture  boundary  into  upper  and  lower 
branches,  and  for  each  branch  into  breakup  rfnd  fill-back 
portions. 


3  Two  Sliding  Surfaces 

There  are  self-acting  bearings  with  two  sliding  surfaces.  For 
instance,  both  the  shaft  and  the  bushing  may  be  ro:a  ing,  each 
at  a  different  speed,  while  their  centers  are  kept  at  a  fixed 
displacement  from  each  other.  The  applicable  equations  for 
such  a  bearing  are 

h.c-, 


:i<: £)($]♦»[(/ snsk--.*--*] 


*0 


( 1| .  I;)  -  (u!) ,  <  22 a.  b,  c) 

In  older  to  describe  the  adhered  film  in  the  cavitation 
domain,  it  is  necessary  to  assume  a  f  rame  of  reference  such 
that  one  of  the  surfaces  does  not  have  a  sliding  velocity.  Such 
a  frame  of  reference  may  be  rotating  at  u>,,  and  in  such  a 
frame  of  reference,  the  bearing  gip  is  no  longer  time- 
independent.  The  required  transformation  is 

x =A-P,/  (23) 

Accordingly,  equations  (22<i.  b)  become 

h  =  C-e  cos(-^-) 


The  difference  between  equations  (22 b)  and  (24 b)  is  only  in 
appearance  so  far  as  the  pressure  field  is  concerned.  However, 
the  Couettc  flux  now  becomes  1/2  {V2  ~  V\)h  so  that  the 
sweeping  motion  of  the  rupture  boundary  is  governed  by 


+  (V3-V,)\ 


(,-*-)] 


sin/* 


(25) 


0,  is  the  sweeping  speed  viewed  in  the  rotating  frame  of 
reference.  tiu  is  the  thickness  of  the  adhered  film  which  is 
attached  to  the  moving  surface  (in  the  rotating  frame  of 
reference).  And.  in  terms  of  the  coordinates  of  the  rupture 
point. 


-  [(£><SP.(S)-(S)(S). 
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(h,  -  h  t  > 


At 


CJO.'IOO. 

"O':  -  *'.>(*■  -A.)  f:  ]  (2<w.M 


Ai  the  veadv -state  condition,  /— »,  the  rupture  boundary 
should  be  .vine  stationary  in  the  space  fixed  coordinate 
system,  i.e  , 


lim 


Mr 

dl 


Accordingly,  equation  (26t/)  becomes 


lim 

;  -  jo 


h,' 

12,. 


IMtHt). 


♦  *  ( +  Vi)h,  -  I  (21) 

At  the  b  cakup  point,  equation  (27)  is  satisfied  by 


Since  h{  cannot  exceed  hh%  it  is  nccessar\  to  impose  the 
constraint 

ir,l>U'1»  (29) 

This  means  the  choice  of  the  rotating  frame  of  reference  is 
not  arbitrary,  it  must  Ik*  fixed  to  the  slower  surface.  Ac¬ 
cordingly,  t he  adhered  film  must  he  attached  to  the  faster 
surface. 

Having  established  tb„*  “birth  value”  of  the  adhcied  film 
and  its  transport  spcec1,  it  is  not  necessary  to  retain  the 
rotating  fiamc  of  reference  in  subsequent  analysts.  One  can 
now  retun  to  the  space -fixed  coordinate  system  to  describe 
the  sweeping  motion  of  t  ie  rupture  boundary: 


sin,ij 

+  [|(>e  I'a/r,  -  y./rjsimi) 

(30) 

where 

3  =  cos  1  ( V  *  *n./  k':)  ±  r /2 

*-<*-«  ((&)[($.(£ 

) 

> 

(31) 

(£>. 

] 

(32) 

At  a  suKiitical  breakup  point 

£).] 

sKj 

(33) 

The  adhered  film  is  fo  med  with  its  “birth”  thickness 

(34) 

At  a  sut  a  weak  up  point,  h  of  prior  time  applies  in 

equations  *  2).  The  characteristic  equation  ia  the 

cavitation  dot  an  is 

v(M  -  \10>=  l\t  (35) 

filiations  (.10  35)  are  seen  to  include  the  “convent  onal" 
self-acting  tearing  (fi,  *  0.  k:  =  H  as  a  special  case, 
corresponding  to  equations  (5,  8,  9.  12.  1  lo.  10).  respectively. 

The  counter-rotating  bearing  ( f,  =  l':)  without  squeeze 

is  of  no  interest.  The  film  pressure  does  not  develop  any 
circumferential  variation,  and  therefore  its  film  docs  not 
rupture. 

The  corotational  bearing  (!’,  =  l\)  vields  the  special  birth 
thickness  of  the  adhered  film 

(36. 

and  the  sweeping  motion  at  a  breakup  point  is  preempted  b> 
the  Sw  ift -Slieber  condition 


4  Squeeze  f  ilm  without  Sliding 

The  pressure  of  a  pure  squeeze-film  satisfies 


and  a  suitable  set  of  external  boundary  conditions,  e.g., 
equation  (20).  Rupture  may  also  take  place  as  a  moehai  isin  u 
limit  the  lower  bound  of  the  squeeze -film  pressure.  Because 
the  surfaces  have  no  sliding  motion,  fluid  left  in  the  ca\  itatior 
domain  is  “trapped.”  The  motion  of  the  rupture  boundary  i* 
given  by  equation  (4)  upon  selling  V  to  zero: 


The  geometry  of  the  cavitation  domain  can  vitl  be 
icpresented  b>  1  ig.  I;  however,  in  the  absence  oi  the  sliding 
velocity,  the  concepts  of  breakup  and  fill -bask  must  be 
reexamined. 

If  the  film  pressure  is  above  the  cavitation  pressure,  then 
the  right-hand  side  is  positive.  Since  the  cavitation  domain  i 
generally  not  filled,  (h,  -/»„)  would  be  positive,  so  t  up  V 
would  also  be  positive.  Thus  the  rupture  boundary  woild 
move  into  the  cavitation  domain  and  the  trapped  fluid  o' 
prior  time  would  become  a  part  of  the  squeeze  film.  Tiis  is  ;t 
fill-back  point,  and  htl  represents  the  amount  of  napped  fluid 
already  existing  on  the  cavitation  side  of  the  rupture  boun 
dary. 

Under  certain  circumstances,  the  constraint  on  the  lowe* 
bound  of  the  film  pressure  may  dictate  that  U,  be  negaiive.  si* 
that  the  rupture  boundary  would  recede,  leaving  behind  new 
fluid  to  become  trapped  in  the  cavitation  domain.  But  the 
right-hand  side  of  equation  (39)  must  be  nonnegative 
Consequently,  h ,  must  equal  the  film  thickness  to  rc:orcile 
the  two  sides  of  equation  (39),  while  the  Swift- Stieber  con 
dilion  is  to  be  satisfied  by  the  squeeze  Him.  The  receding 
rupture  boundary  of  a  squeeze  film  can  therefore  be  regarded 
as  a  breakup  point.  The  corresponding  birth  value  of  is 

hy  =  tih\  [Vp)h~ 0  if  Vh< 0  (RV /)i 

It  is  possible  for  a  pure  squeeze-film  to  have  a  fully  receding 
rupture  boundary,  rhis  usually  happens  when  the  film 
thickness  is  increasing  with  time  for  the  entire  cavitaton 
domain.  Thus  one  can  have  a  filled  cavitation  domain  the 
initial  state,  so  tong  as  the  Swift-Stiebcr  condition  is  enforced 
on  the  entire  rupture  boundary  in  the  film  pressure . 

The  half  Sommer  fold  film  pressure  may  also  be  used  as  the 
initial  condition.  When  this  is  done,  a  partially  scavenged 
initial  cavitation  domain,  h (r,  c,  0)  <  h  (*,  2,  0).  must  be 
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assumed.  Otherwise  U,  would  be  infinite,  rendering  the  initial 
state  meaningless. 

To  consider  a  whirling  squeeze-film,  equation  (39)  can  be 
replaced  by  the  components 


-(g)[(D. -(*),(*)] 


(4la.6) 

Upon  reaching  a  steady-whirl  condition, 

i  dx'  i- 

lint  — —  =  l 

I  --  jo  dt 

then  equation  (41  a)  becomes 

which  is  satisfied  by 


IS  (%).•( 


(v 


dp  > 


lim 


(  s' )  ‘ :- 


(I), 


(42*.  b) 


These  relations  represent  a  precise  analog  of  the  steadily 
whirling  squeeze-film  to  a  statically  displaced  selfacting 
journal  in  both  the  film  pressure  and  the  rupture  boundary. 
However,  in  the  latter  case,  hu  -  1/2  hh%  i.e.,  the  fluid 
content  in  he  cavitation  domain  of  the  journal  bearing  is  one- 
half  of  tha:  of  the  whirling  damper. 


5  Subcavitation  Film  Pressure 

The  assumption  of  a  constant  and  uniform  pressure  in  the 
cavitation  domain,  p<.  does  not  necessarily  mean  it  has  to  be  a 
lower  bound  of  the  film  pressure  at  the  same  time.  The  latter 
condition  should  be  recognized  as  a  conceptual  artifact  which 
is  invoked  to  avoid  a  rigorous  consideration  of  the  ther¬ 
modynamic  and  chemical  processes  which  may  be  present. 
Such  a  posture  can  be  justified  if  one  is  concerned  with  a 
steady-state  situation  and  that  thermodynamic  and  chemical 
equilibria  tan  be  assured.  Unfortunately,  thermodynamic  and 
chemical  equilibria  hardly  ever  exist  in  the  pure  sense.  In 
studying  dynamic  problems,  it  is  quite  possible  that  the  rate  of 
chemical  k  netics  associated  with  mass  transfer  between  liquid 
and  gaseous  phases  in  a  rupturing  film  may  be  quite  small  in 
comparison  with  the  rate  of  mechanical  phenomena  of  in¬ 
terest.  Therefore,  it  ma>  be  necessary  to  abandon  the  con¬ 
straint  that  p  is  the  lowet  bound  of  the  film  pressure. 

One  approach  to  resolve  the  conceptual  dilemma  described 
earlier  is  to  specify  the  existence  of  a  cavitation  domain,  in 
which  p€  prevails,  but  not  to  require  a  lower  bound  of  the  film 
pressure  in  the  unruptured  domain.  Such  a  model  may  be 
justified  by  the  following  lines  of  reasoning. 

•  Void  generation  in  the  fluid  film  begin  at  internal 
nuclcation  sites  of  entrained  bubbles,  as  well  as  surface 
nucleation  sites. 

•  There  are  sufficient  number  of  nucleation  sites  where  ihe 
growth  rate  is  high  enough  so  that  the  voids  merge,  aggregate, 
and  cumulate  within  the  boundaries  of  the  bearing  film. 

•  After  an  adequate  lapse  time,  which  is  short  in  com¬ 
parison  with  the  operation  time  of  the  bearing  but  long  as 
compared  with  the  growt  It -aggregation-cumulation  processes, 

• 


mass  transfer  between  the  phases  become  minimal  and  can  be 
neglected. 

•  In  the  absence  of  significant  mass  transfer,  rcsidu il  void 
fraction  in  the  fluid  film  is  negligible.  Pressure  fluctuations 
within  the  fluid  film  can  be  sustained.  In  any  ease.  with  the 
presence  of  a  cavitation  domain,  the  pressure  within  the  fluid 
film  cannot  Tall  substantially  below  pK . 


On  the  basis  of  these  postulated  ideas,  one  can  proceed  to 
examine  the  mathematical  consequences  of  subcavitation 
pressure  in  the  fluid  film  near  the  rupture  boundary. 

Consider  first  the  self-acting  bearing.  With  a  subcavitation 
film  pressure  in  the  immediate  vicinity  of  the  rupture 
boundary,  ( V/?)r*n,  would  be  positive  and  the  Poiseuille  flu\ 
would  aim  away  from  the  cavitation  domain.  The  total  flux 
toward  a  rupture  boundary  is 

<?.=  [-(^)(V/j)r+^  (V,  +Vj)A,].«,  (431 

r  may  be  replaced  by  b  or/ in  reference  to  a  breakup  or  a  fill 
back  point.  With  a  subcavitation  film  at  a  breakup  point,  Q, 
can  be  negative.  When  this  happens,  the  breakup  boundary 
withdraws,  leaving  behind  a  dry  region.  Therefore,  a  drying 
up  condition  of  a  breakup  point  is  due  to 

<?/><  0 

consequently . 

hu-  0;  Uh=Qb/hh<  0  (44</, />) 

If  Qh  is  positive,  but  less  than  the  Couette  flux  due  to  a 
subcav  itation  film,  then  the  adhered  film  thickness  .x  si  ill 
given  by  equation  (34);  in  the  meantime,  as  may  be  obtained 
from  equation  (30), 

<J5> 

indicating  a  receding  motion.  If  the  film  pressure  at  a  flll-back 
point  is  subcavitation,  then  a  supercritical  fill-back  condition 
can  exist: 


((*>Ks),(£)-(S)J 


+  [  *  (  V,  +  »',)/),  -  )  =  r.  +  ±1  >  I  (46) 

The  fill-back  point  would  be  moving  ahead  the  transported 
adhered  film,  opening  up  a  dry  domain;  accordingly,  the 
supercritical  sweeping  speed  of  a  fill -back  poin;  is  obtained  by 
setting  to  zero  in  equation  (30)  yielding 

U,  =  Q,/h,< 0;  where  Q,fh,  ..sintfjV,  +  (  ).u  j 

(47) 


h9  being  given  in  the  prior  time  value. 

The  pure  squeeze-film  with  subcavitation  film  pressure  can 
be  regarded  as  a  special  case  of  the  self-acting  bearing  w  th  the 
sliding  velocities  set  to  zero.  The  total  flux  is  then  identically 
Ihe  Poiseuille  flux, 

/  h  *  \ 

(?,  =  -(i~)<Vp)r.n,  (48) 

and  there  is  no  distinction  between  breakup  and  fill-back 
points.  The  subcavitation  film  pressure  automatically  rauscs 
the  drying  up  condition: 

ha  =0;  Ur  =  Q,lhr  if  Q,<0  (49u,  b) 

The  possibility  of  dry  regions  due  to  the  subcavitation  film 
pressure  may  be  described  in  the  world  diagram  as  shown  in 
Fig.  4.  The  world  diagram  of  a  pure  squeeze-film  has  vertical 
characteristic  lines. 
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Fig.  4  World  diagram  with  aubcavHatlon  film  pressure 


6  In  Retrospect 

The  lime  dependent  problem  of  film  rupture  can  be 
analyzed  in  terms  of  the  movement  of  the  rupture  boundary. 
The  instantaneous  motion  of  the  rupture  boundary  is  the 
result  ol  Tow  balance  between  the  full  film  on  one  side,  which 
is  described  by  the  traditional  lubrication  theory,  and  the 
ruptured  film  on  the  other  side,  which  is  incapable  of 
sustaining  any  significant  pressure  gradient  due  to  the 
abundant  presence  of  large  scale  voids. 

Consistent  with  the  lubrication  analysis  employed  for  the 
full  film,  inertia  effects  may  be  neglected  for  the  ruptured 
film.  Using  a  homogeneous  model,  the  ruptured  film  may  be 
regarded  as  a  layer  adhered  to  the  sliding  surface.  The 
adhered  film  is  transported  through  the  cavitation  domain 
with  an  nvariant  thickness  which  cannot  exceed  the  local  gap. 

The  adhered  film  next  to  a  ruptured  boundary  may  have 
one  of  three  states,  in  reference  to  its  thickness.  The  adhered 
film  is  in  its  prior  tin  e  state ,  if  the  rupture  boundary  is  ad¬ 
vancing  on  the  sliding  surface;  the  thickness  of  the  prior  time 
adhered  film  is  determined  from  past  history  of  the  ruptured 
film.  The  adhered  film  is  in  its  newly  born  state  if  the  rupture 
boundary  is  "subcritually  fed”  by  the  full  film  flux  and  lags 
behind  i  he  sliding  sutfacc;  the  thickness  “ui  birth”  of  the 
adhered  film  is  one-half  of  the  local  gap  (if  only  one  of  the 
two  surfaces  has  a  sliding  speed).  The  adhered  film  is  in  its  dry 
state  if  i  he  total  film  flux  is  directed  away  from  the  rupture 
boundary;  the  dry  adhered  film  has  no  thickness.  The  dry 
adhered  film  always  borders  a  full-film  with  a  subcavitation 
pressure. 

In  mathematical  teims,  dynamics  of  the  ruptured  film  is 
concerned  with  a  mixed  field  initial  value  problem.  In  the 
model  problem  illustrated  in  this  paper,  the  full-film  domain 
is  the  outer  field;  the  field  variable  is  the  film  pressure,  its 
governing  equation  is  elliptical  in  character.  The  cavitation 
domain  is  the  inner  field;  the  field  variable  is  the  adhered  film 
thickness,  its  governing  equation  is  hyperbolic  in  character. 
Coupling  between  the  two  fields  is  maintained  through 

•  the  cavitation  pressure  as  the  internal  boundary  con¬ 
dition  of  the  film  pressure, 

•  birth  of  new  adhered  film  at  subcritically  fed  rupture 
point,  and 

•  motion  of  the  rupture  boundary  as  determined  by  the 
states  of  the  fields  on  both  sides. 

The  initial  value  set  includes  the  location  of  the  rupture 
boundary  and  the  distribution  of  the  adhered  film.  The  film 
pressure  field  at  every  instant  can  be  computed  according  to 
Reynolds'  equation  with  the  cavitation  pressure  imposed  at 


the  rupture  boundary.  Birth  values  of  adhered  film  thickm> 
as  well  as  the  sweeping  speed  of  the  rupture  boundary,  can  t  c 
determined  accordingly.  Time  domain  integration  of  ve 
motion  of  the  rupture  boundary  allows  incremental  update  g 
of  both  field  variables  and  the  rupture  boundary.  II  the  fil  n 
pressure  is  not  allowed  to  become  lower  than  the  cavitatii  n 
pressure,  the  initial  state  of  the  adhered  film  cannot  he  a  ful  > 
filled  gap. 

Generalization  to  the  case  of  two  sliding  surfaces  requir  % 
that  the  ruptured  film  be  adhered  to  the  faster-sliding  surfa*.  \ 
The  birth  value  of  the  adhered  film  is  thicker  b>  an  addition  il 
fraction,  which  is  equal  to  the  ratio  of  the  slower  sliding  spe  d 
over  the  faster  sliding  speed,  as  compared  to  that  for  the  ca  e 
of  a  single  sliding  surface.  The  adhered  film  is  transjioi  ted  1  \ 
the  faster  sliding  speed.  The  case  of  equal  sliding  speed,  in¬ 
cluding  the  pure  squeeze  film  with  zero  sliding  of  both  su 
faces,  has  a  birth  value  of  the  adhered  film  equal  to  the  ti  II 
local  gap,  and  the  Swift-Stieber  condition  is  to  be  satisfied  a>  a 
subcritically  fed  rupture  point. 
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